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ABSTRACT 

We present the rest-frame near-ultraviolet (NUV) spectroscopy of star-forming galaxies (SFGs) at 
0.6 < z < 1.2 from the Extended Baryon Oscillation Spectroscopic Survey (eBOSS) in SDSS-IV. One 
of the eBOSS programs is to obtain 2" (about 15 kpc) fiber spectra of about 200,000 emission-line 
galaxies (ELGs) at redshift z > 0.6. We use the data from the pilot observations of this program, 
including 8620 spectra of SFGs at 0.6 < z < 1.2. The median composite spectra of these SFGs at 
2200 A < A < 4000 A feature asymmetric, preferentially blueshifted non-resonant emission, Fe II*, 
and blueshifted resonant absorption, e.g., Fe II and Mg II, indicating ubiquitous outflows driven by 
star formation at these redshifts. For the absorption lines, we find a variety of velocity profiles with 
different degrees of blueshift. Comparing our new observations with the literature, we do not ob¬ 
serve the non-resonant emission in the small-aperture (< 40 pc) spectra of local star-forming regions 
with the Hubble Space Telescope (HST ), and find the observed line ratios in the SFG spectra to be 
different from those in the spectra of local star-forming regions, as well as those of quasar absorption¬ 
line systems in the same redshift range. We introduce an outflow model that can simultaneously 
explain the multiple observed properties and suggest that the variety of absorption velocity profiles 
and the line ratio differences are caused by scattered fluorescent emission filling in on top of the ab¬ 
sorption in the large-aperture eBOSS spectra. We develop an observation-driven, mo del-independent 
method to correct the emission-infill to reveal the true absorption profiles. Finally, we show that the 
strengths of both the non-resonant emission and the emission-corrected resonant absorption increase 
with [O II] AA3727, 3730 rest equivalent width and luminosity, with a slightly larger dependence on 
the former. Our results show that eBOSS and future dark-energy surveys (e.g., DESI and PFS) will 
provide rich datasets of NUV spectroscopy for astrophysical applications. 

Subject headings: surveys - ultraviolet: galaxies - galaxies: evolution intergalactic medium - galaxies: 

ISM - galaxies: star formation - quasars: absorption lines 


1. INTRODUCTION 

Spectroscopy is one of the most important tools in 
modern astronomy. Historically, ground-based observa¬ 
tions have mostly been focused on the optical window 
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of the spectrum because of the extinction effects of the 
Earth’s atmosphere and the low brightness of most cos¬ 
mological sources. In the near-ultraviolet (NUV) range 


at 2000 A < A < 3000 A, the most-studied sources 
are quasars (e.g., iSchmidtl 119631 : iFrancis et all 119911 : 
iVanden Berk et al.ll200lD and quasar ab s orption-line sys¬ 
tems (e.g., iBurbidge fc Kinmanl 119661 iBergeronl 119861 : 


iSteidel fe Sargentlll992t iChurchill et alil2000( ) at moder¬ 

ate redshifts (0.5 < z < 2.5), thanks to the high energy- 
conversion efficiency of supermassive black holes. Af¬ 
ter space-based telescopes became available, such as the 
International Ultraviolet Explorer ( IUE) and the Hub¬ 
ble Space Telescope (HST), astronomers started to ex¬ 
plore the UV part of the spectral energy distributions 
(SEDs) of stars and galaxies. Most of the space-based 
observations are dedicated to the far-UV (FUV) window 


at A < 2000 A because of the large number of transi¬ 
tions available at those energies, leaving the NUV win¬ 
dow of the SEDs of stars and galaxies largely-unexplored 
territory. Most of our knowledge about galaxy SEDs 
in the NUV have originated from most recent observa¬ 
tions with the 8-meter or large r telescopes such as the 


Keck and VLT telescopes fe.g.. iMartin & Bouche 

2009 

Weiner et al. 2009: Rubin et al. 20lIF Talia et ah 

2012 

Erb et al.1120121). 


The most useful spectral features in the NUV 
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have been the resonantPl absorption lines induced 
by the low-ionization species Mg II and Fe II in 
the interstellar medium (ISM) and/or circumgalac- 
tic medium(CGM). In the spectra of star-forming 
galaxies (SFGs), these absorption lines are observed 
to be almost always blueshifted from the systemic 
velocity of the galaxy, which serves as direct ev¬ 
idence f or ubiquitous outflows driven by star for- 


rnation (Tremonti et a L| 20071 

Martin & Bouche 2009t 

Weiner et al. 200Sf Rubin et al. 

201(1 2014 

; Talia et al. 

2012: Erb et al. 20ll Martin et al. 2012: 

Kornei et al. 

2011)- 


More recently, observations have also revealed the non¬ 
resonant emission features Fe II* in the NUV spectra 
of SF Gs dRubin et all 1201 ll iTalia et ahll2012t lErb et all 
IMl IKornei et all 120121 I201.C These non-resonant 
emission lines are likely caused by scattered fluores¬ 
cent photons after the occurrence of resonant absorp¬ 
tion. As the non-resonant emission lines i n the FUV 
(e.g., iShaplev et al.l[200l iJones et al.l 120121) . they have 
provided a new means to study the radiative transfer 
physics in the complex baryon processes in galaxy evo¬ 
lution (e.g. IRu bin et al.1 1201 It IProchaska et all 120111 : 
IScarlata fe Panagiall2015f) . 

The non-resonant emission, e.g., Fe II*, was sometimes 
referred to as fine-structure emission in the literature, 
because the lower energy level of the transition is an ex¬ 
cited level due to fine-structure splitting, as opposed to 
the ground level for the corresponding resonant transi¬ 
tion. The term fine-structure emission , however, is used 
more often for the emission between two fine-structure 
levels (with different J-values) in the same term (with 
the same L , S), which usually occurs in the infrared (IR), 
such as [C II] A157.7 pm. To avoid confusion, we here 
simply use the term non-resonant emission and reserve 
fine-structure emission for those infrared emission lines. 
We also label the (permitted) non-resonant emission lines 
with a right superscript asterisk “*”, although the ori¬ 
gin of this convention appears to have been lost in the 
literaturcQ. 

The fine-structure splitting causes multiple transitions, 
i.e., multiplets, with energies only slightly different from 
each other . In her pioneering multiplet tables dMoorel 
119451 [T950h . Charlotte E. Moore grouped the multiplets 
in the ascending order of the excited terms and this num¬ 
bering convention has been widely used since. We follow 
t he same co nvention for the multiplets in the UV ( Moore 
I1950L 1952 b though we skip the numbering of multiplets 
in the optical as in the recent literature. In Appendix [Al 
we summarize our choices of conventions and also present 
the atomic data needed in this analysis for references. 

In the last decade, the detection of the baryon acoustic 
osci llation (BAO) in the la r ge-scale distribut ion of galax¬ 
ies (lEisenstein et al.l 120051 [Cole et all 120051 ) has shifted 
the focus of observational cosmology to massive redshift 
surveys over the entire observable sky, and we are now in 

13 We refer to a transition between two terms as a resonant 
transition if and only if the lower energy level is the lowest (ground) 
level of the ground term. See Appendix El for clarification. 

14 In SYMBOLS, UNITS, NOMENCLATURE AND FUNDA¬ 
MENTAL CONSTANTS IN PHYSICS, or the “Red Book”, com¬ 
piled by International Union of Pure and Applied Physics (IUPAP), 
the authors reco mmended this conv ention for excited atomic levels 
(see Section 2.1. 1 Cohen et ahll 19871) . 


the course of measuring the scale of BAO as a function of 
cosmic time and chronicling the expansion history of the 
Universe. Among them are the Exte nded Baryon Oscilla- 
tion Spectroscopic Survey (eBOSS. iDawson et al.ll2015l) 
within the fourth phase of the Sloan Digital Sky Survey 
(SDSS-IV, Blanton et al. in prep.), the D ark Energy 
Spectroscopic Instrum ent survey (DESI, iSchleeel et al.l 
1201 It iLevi et all I2013D. and the Prim e Focus Spectro¬ 
graph survey (PFS. lTakada et akll2014D . These programs 
will obtain optical spectra in the observer frame for tens 
of millions of galaxies at redshift 0.6 < 2 : < 2.0, provid¬ 
ing unprecedentedly rich spectroscopic datasets with the 
coverage in the rest-frame NUV. 

Among the new BAO projects, the eBOSS survey 
started taking data in Fall 2014. One of the large- 
scale structure tracers in eBOSS is emission-line galaxies 
(ELGs). We have obtained about 12,000 ELG spectra 
across the redshift range 0 < 2 : < 1.5 in the pilot obser¬ 
vations, conducted to optimize targeting strategies. We 
discuss the whole sample briefly and focus on the NUV 
spectroscopy of 8620 ELGs at redshift 0.6 < z < 1.2. 
We describe the observations and the dataset in detail in 
Section [2] 

The rest of the paper is organized as follows. In Sec¬ 
tion^ we present the composite spectrum of ELGs in the 
NUV and compare it with those of quasar absorption-line 
systems and local star-forming (SF) regions. We inter¬ 
pret the observations with a spherical gas outflow model 
in Section [1] In Section [5] we investigate the dependence 
of the spectral features on the [O II] AA3727, 3730 proper¬ 
ties. We summarize our results in Section [H] When nec¬ 
essary, we assume the ACDM cosmogony, with Ua = 0.7, 
= 0.3, and Ho = 70kms -1 Mpc -1 . Throughout the 
paper, we use vacuum wavelength and some of the lines 
are therefore labeled as 1 A longer than if labeled with 
the air wavelength, such as [O III] A5008. 

2. DATA 

2.1. The eBOSS Survey 

The eBOSS survey is the cosmology survey 
within SDSS-IV, an extension of th e BOSS survey 
(IDawson et alJ 1201311 in SDSS-III (lEisenstein et all 
(Ml . Over a six-year period beginning in Fall 2014, 
eBOSS will observe four independent tracers of the 
underlying density field: luminous red galaxies (LRGs), 
quasars, Lya forest, and ELGs, and measure the cosmo¬ 
logical distances as a function of redshift with BAO as a 
standard ruler and record the expansion history of the 
Universe at 0.6 < 2 : < 2.3. 

The eBOSS survey uses the sam e two i dentic al, multi¬ 
object spectrographs as in BOSS dSmee et all 12013T) on 
the 2.5-meter SDSS Telescope (|Gunn et al.ll2006[ ) at the 
Apache Point Observatory in New Mexico. Each spec¬ 
trograph is equipped with two cameras, one blue and one 
red, with a dichroic splitting the light at roughly 6000 A 
and a combined coverage spanning between 3600 A and 
10400 A. The spectral resolution 1Z is 1560 — 2270 in the 
blue and 1850 — 2650 in the red channel, with a mean 
approximately 7 Z « 2000. The spectrographs are fed by 
1000 optical fibers (500 for each), covering a field-of-view 
(FoV) of about 7.5 square degrees. The aperture diame¬ 
ter of the fibers is 2 ", smaller than 3 " in SDSS-I/II. The 
typical total exposure time for each pointing is about 75 
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minutes. 

eBOSS selects targets primarily based on the 
SDSS imaging, obtained through a set of ugriz fil¬ 
ters (iFuk ugita et all I199 (tT) with a wide-held camera 
(IGunn et al.1 llOOSF in a drift-scan mode. As the de¬ 
sired targets are fainter than those in SDSS I-III, in 
order to achieve high targeting efficiency, eBOSS in¬ 
cludes supplementary imaging data with other_ instru¬ 
ments, including the infrared photometry (lLang et al.l 
2014fl bv the Wide- Held Infrared Survey Explorer ( WISE , 
Wright et al.ll201(l f . the U-band imaging from the South 
Galactic Cap U-band Sky Survey (SCUSSJ[3 conducted 
at the 2.3-meter Bok telescope at the Kitt Peak Na¬ 
tional Observatory, and the de ep qrz imaging with the 
Dark Energy Camera (DECam, iFlaugher et al.ll2012f) on 
the 4-meter Blanco telescope at the Cerro Tololo Inter- 
American Observatory (CTIO) in Chile. 

The spectral-reduction and redshift-htting pipeline 
in eBO SS is a c ontinuation of the BOSS pipeline 
(|Bolton et al.l l2012h improved to yield required perfor¬ 
mance on spectra with lower signal-to-noise ratio (S/N) 
than in BOSS. As BOSS, eBOSS classifies objects and de¬ 
rives redshifts based upon principal-component-analysis 
(PCA) templates for quasars and galaxies, and archety¬ 
pal templates for stars, though the team is currently ex¬ 
ploring an alternative approach fully based on archetypal 
templates. For the analysis presented here, we use the re¬ 
sults derived from the pipeline version v5_7_8, which still 
uses the PCA templates for galaxies and quasars. The 
pipeline flags a redshift with a warning (ZWARNING) based 
on various quantitative criteria, such as a small \ 2 differ¬ 
ence between the best and second-best fits. Tests have 
shown that the ZWARNING is conservatively defined and 
the redshift success rate for objects with ZWARNING==0 
is better than 99%. We only use objects classified as a 
GALAXY with ZWARNING==0 in this analysis. 

For more details regar ding the eBOSS survey, we refer 
the reader to lDawson et alJ (l2015h . 

2.2. Emission-line Galaxies 

Emission-line galaxies are one of the four tracers 
targeted in eBOSS. The primary tracer of the large- 
scale structu res in BAO observatio ns (at z < 1) has 
been LRGs (lEisenstein et al.l 120011) . because of their 
broadband brightness and well-understood SEDs. At 
higher redshift, however, most of the distinctive fea¬ 
tures of LRGs, such as the 4000 A break, G4300- 
band and MgH/Mg b absorption bands, are redshifted 
into the Meinel hydroxyl forest. Blueward of the 
4000 A break, the cool giant stars dominating the con¬ 
tinuum emit little flux, which makes spectral classi¬ 
fication and redshift determination difficult. More¬ 
over, the number density of red galaxies is l ower at 
highe r redshift due to ga l axy evolution (e.g . . iBell et al.1 
12004 iFaber et~ahl 120071 : IMoustakas et al.1 I2013IL ren- 
dering them less useful for cosmological purposes. 
On the other hand, the cosmic star-formation rate 
(SFR) density increases precipitously with redshift, and 
at z ~ 1, is about a n order-of-magnitu d e higher 

than at present (e.g., Madau et al.1 Il996t iHopkinsI 
12004 IMadau fe Dickinsonl 201411 . Star-forming galaxies 

15 http://batc.bao.ac.cn/Uband/ 


(SFGs) exhibit the strong [O II] AA3727,3730 emission 
feature, which is a doublet distinguishable at a medium 
spectral resolution 7 Z ~ 4000. The number density of 
[O II] AA3727, 3730 emi tters increases steadily with red- 
shift up to z ~ 2 (e.g.. IZhu et a.l.l[20091 IComparat et al.l 
120154 ISobral et al. I2015lf. The next-generation BAO 
surveys, e .g., DESI (Schlegel et al.l20lHlLevi et al.l[2013h 
and PFS dTakada et al.ll20l41 . will primarily target these 
[O II] emitters at redshift z > 1. 

At redshift z < 1, the strong emission lines of SFGs 
enhance their brightness in the optical (in the observer 
frame), facilitating their detection even in the r elatively 
shallow SDSS imaging. IComparat et al.1 (|2013H demon¬ 
strated that it is feasible to conduct an ELG survey with 
the BOSS spectrographs that can achieve the number 
density and the volume coverage required for BAO mea¬ 
surements at 0.6 < z < 1.0, with targets selected from 
the SDSS imaging. As the continuum of SFGs is domi¬ 
nated by the emission of hot O/B stars in the rest-frame 
NUV and that of dust and polycyclic aromatic hydro¬ 
carbons (PAHs) in the IR, u/U -band and/or IR data 
can help improve the targeting efficiency of these ob¬ 
jects. We have been working on optimizing the ELG se¬ 
lection strategies with additional data, including the U- 
band photometry from SCUSS, the IR photometry from 
WISE , and the deeper grz imaging with DECam. 

Favorable weather during SDSS-III led to an early com¬ 
pletion of the BOSS survey. A fraction of the remaining 
time was allocated to an eBOSS pilot program known 
as the Sloan E xtended Quasar, ELG, and LRG Sur¬ 
vey 1SEQUELS. lAlam et ahlfeOlSf) . In Fall 2014, within 
eBOSS, we also conducted a series of pilot observations 
to test possible techniques for the ELG target selection. 
With data from these pilot observations, the team is 
currently investigating different selection algor ithms to 
maximize the targeting efficiency dRaichoor et al.l 120151: 


argetmg efficiency (Kaiciroo' ct al. ZUIo; 

1l2015btlDelubac et al.ll2015HJouvel et al.1 


Comparat et al 


2015T) . The ELG cosmology survey will begin in Fall 
2016, the third year of eBOSS, with the optimal selection 
strategy to be defined from these pilot data and investi¬ 
gations. The survey aims at obtaining secure redshift for 
about 200, 000 ELGs at 0.6 < z < 1.0 and measuring the 
BAO scale with an accuracy of about 2% at an effective 
redshift (z) ~ 0.8. 

For more details regarding the ELG target selec¬ 
tions and the cosmological appli cations, w e re fer the 
reader to reference s above and lDawson et al.1 (|2015t) and 
iZhao et al. I (|2015j h 


2.3. The ELGs from the eBOSS Pilot Observations 

In total, the BOSS/SEQUELS ancillary program and 
the early eBOSS pilot observations provided about 
12,000 ELGs spanning 0 < z < 1.5, peaked at z ~ 0.8. 
We show the redshift distribution in Figure [T] In Fig¬ 
ure El we present the median composite spectrum of all 
the ELGs, in the wavelength range 2000 < A < 7500 A. 
As we included all the ELGs in this composite spec¬ 
trum, the objects contributing at each wavelength are 
different, with high-redshift objects dominating in the 
NUV and low-redshift ones in the optical. The compos¬ 
ite spectrum appears as a typical active SFG spectrum 
(e.g., lKennicuttlll992| h We have labeled the prominent 
features in the figure. In the optical, it features strong 
hydrogen Balmer recombination emission on top of rel- 
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Rcdshiff Distribution 



Z 

Fig. 1.— The redshift distribution of all emission-line galaxies 
(ELGs) from the eBOSS pilot observations. The solid area shows 
8620 galaxies at 0.6 < z < 1.2. 

atively weak stellar Balmer absorption, strong nebular 
forbidden lines due to collisionally-excited metal atoms, 
e.g., N, O, N+, 0+, S+, 0 2+ , Ne 2+ , and Ar 2+ . As the 
BOSS spectrographs cover the [O III] A5008 up to red- 
shift ^ ~ 1.0, we expect that eBOSS will provide an 
important opportunity for studies of the ISM properties, 
such as the gas-phase metallicity, in strong SFGs at mod¬ 
erate redshift. The full composite spectrum also displays 
some relatively weak stellar metal absorption lines, e.g., 
the Fraunhofer Ca II H & K, G4300-band, Mg I b, and 
Na I D lines. We summarize the identified lines in the 
list given in Appendix [A] 

Comparat et al. (2015b) describes the full sample in 
detail and we refer the reader to that paper for more in¬ 
formation. The mean precision of the redshift at z ~ 0.8 
is about 20kms -1 thanks to the strength of the emis¬ 
sion lines. From the SED fitting of the composite spec¬ 
tra with different line strengths, the average stellar mass 
(M*) at z ~ 0.8 is about 10 lo M Q . The investigations 
of the physical properties of individual galaxies, such as 
stellar mass (M*), SFR, and metallicity, are ongoing and 
will be presented in future papers. 

The redshift coverage of the ELG observations at z > 
0.6 allows us to probe the NUV part of their SEDs, as 
shown in Figure [2] The composite spectrum between 
2900 A and the Balmer Break at 3646 A is essentially 
featureless except for the weak He IA3189 emission line. 
Blueward of 2900 A, we see strong absorption lines sim¬ 
ilar to those in intervening quasar absorption-line sys¬ 
tems, Mg I, Mg II and Fe II, and non-resonant Fe II* 
emission as well as the C II] emission. The underlying 
stellar continuum rises towards higher energy, as typ¬ 
ical of hot O/B stellar spectra. We fit a power law, 
/(A) oc A^, to the continuum between 2000 A and 2200 A 
and obtain a slope /? ~ —2.1, as expected for an SED 
dominated by O/B stellar emission in the UV. In the 
rest of the paper, we will focus on the absorption and 
emission features in the NUV. 

3. NEAR-ULTRAVIOLET SPECTROSCOPY 


3.1. The NUV ELG Sample 

The eBOSS pilot observations have obtained spectra 
for a large sample of ELGs at redshift z > 0.6 with rest- 
frame NUV coverage, providing us with a good opportu¬ 
nity to investigate the gas processes associated with these 
objects. We here focus on the wavelength range between 
2200 A and 4000 A. We choose the shorter wavelength 
limit so as to cover Fe IIA2250 and Fe IIA2261, and the 
longer limit to cover [O II] AA3727,3730. To ensure the 
same wavelength coverage for all the objects and thus 
the same contributing objects at all the wavelengths, we 
select ELGs between redshift 0.6 and 1.2, which include 
8620 objects. 

As discussed in Section EU we use the reduction out¬ 
puts based on the spectroscopic pipeline version v5_7_8, 
and consider only objects classified as a GALAXY with 
no redshift warning, i.e., ZWARNING==0. The classifi¬ 
cation selection automatically rejects broad-line active 
galactic nuclei (AGNs). We do not make further cuts 
in our sample selection. At the redshifts we are inter¬ 
ested in, some of the lines requir ed in the narrow-lin e 
AGN classification schemes (e.g., iBaldwin et al.lfT981f) . 
such as Ha and [N II] A6584, are not covered by eBOSS. 
For those with [O III] A5008 and H/3 measurements (at 
z < 1) , based on the b lue optical color an d line ratios 
(e.g., lYan et all 120 1 ll [Trouble et all 1201 111 . we expect 
the fraction of narrow-line AGNs to be at most a few 
percent. We also do not expect many lo w-ioniz ation nu¬ 
clear emission-line regions (LINERs. IHecknranlll98C)i ) in 
our sample. The ELGs are selected to be blue galaxies, 
while the majority of LINERs are found in red galaxies 
(e.g., iHo et al.lll997f) . We therefore expect the line emis¬ 
sion in the integrated spectra of the eBOSS ELGs to be 
dominated by contributions from SF activities and will 
use the terms ELGs and SFGs interchangeably. 

3.2. The Method 

The average S/N per pixel in the continuum region 
of the individual ELG spectra is low (< 1) and does 
not allow precise measurements of the absorption fea¬ 
tures for single objects. To study the gas associated with 
the ELGs, we construct high-S/N composite continuum- 
normalized spectra. We use a median estimator, which 
is less prone to extreme outliers. However, we also tested 
our analysis with the arithmetic mean estimator and 
found consistent results, with differences in relative de¬ 
pendences smaller than ler. 

For each observed spectrum, F( A), we first blueshift it 
back to its rest frame on a common wavelength grid. 
We choose the common wavelength grid to have the 
same logarithmic (or equivalently, velocity) spacing as 
in the observer frame, i.e., with dlog 10 A = 1CU 4 or 
du = 69kms _ . In the blueshifting process, we inter¬ 
polate the spectrum with the cubic-B spline method 1 ^!, 
as in the standard SDSS pipeline. We then mask out 
absorption and emission features and fit a cubic polyno¬ 
mial function through the rest of the spectrum. Using 
the best-fit polynomial function as an estimate of the un¬ 
derlying continuum, F’ cont (A), we normalize the observed 
spectrum to obtain the continuum-normalized spectrum 

16 As the spacing is identical before and after interpolation, lin¬ 
ear interpolation yields almost the same results. 


























(F(A)) [arbitrary unit] 



Fig. 2.— The median composite spectrum of all ELGs (at 0 < 2 : < 1.5) at 2000 A < A < 7500 A from the eBOSS pilot observations. We label emission features in green, stellar 
absorption features in brown, and ISM/CGM absorption lines in black. The Greek symbols indicate hydrogen Balmer lines, which appear as nebular recombination emission on top 
of stellar absorption. 
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suited for absorption studies: 


R(X) = 


-Fcont(A) 


( 1 ) 


For a given sample, we construct a median composite 
spectrum of all the continuum-normalized spectra (in 
the rest frame). Finally, we fit a quadratic polynomial 
function to the composite spectrum, again with absorp¬ 
tion and emission features masked out, to remove any 
large-scale residuals, though skipping this final step has 
a negligible effect on the results. We designate the final 
composite as (R( A)): 


{RW) =(£!))• <2) 

Since we mostly work with composite spectra, we will 
drop the ensemble symbol () in the text for simplicity. 

We quantify the absorption and emission strength in 
the continuum-normalized spectra with the rest equiv¬ 
alent width Wo- We define the rest equivalent widths 
for absorption and emission in such a way that they are 
both positive. For absorption, the rest equivalent width 
is given by 

Mmax 

^absorption ^ / [1 - i?(A)] dA , (3) 

•^Amin 


At 2200 A < A < 2900 A, we are able to identify only 
one photospheric absorption line, C III at 2297.58 A. 
However, in subsampling exercises with bootstrapping 
and jackknife, we notice that there are some weak but 
persistent features that are not due to noise (also see 
Section 13.51 below). We do not identify these weak stel¬ 
lar features because the NUV part of the O/B star 
SEDs has not been sufficiently explored in either the¬ 
ory or o bservation. The most recent work was bv 
iRodriguez-Merino et al.1 (120051 the UVBLUE library f^l. 
who built a suite of stellar spectral templates in the NUV 
based on the atmospheric model code ATLAS9 (iKurucd 
119921 : Castelli et al.lll997l l and, for O/B stars, compared 
the model spectra with a few low-resolution spectra taken 
by IUE in 1980s. Although the shape of the underly- 
ing stella r con tinuum is well-understood {j3 ~ —2.0, e.g., 
iKinnev et al.lll993l and Section HOI) , the absorption fea¬ 
tures in theoretical calculations and observations do not 
match each other. For example, we do not detect the 
O III line at 2496 A predicted by the models in our com¬ 
posite sgectrunijjiordid in the IUE observations of O/B 
stars (|Fanelli et al.lll992t IRodrfguez-Merino et al.ll2005l) . 

At 3000 < A < 4000 A, the prominent stellar ab¬ 
sorption features are the well-known Balmer series at 
A > 3646 A, and Ca II H (3969.59 A) & K (3934.77 A) 
lines. The Ca II doublet, however, likely has a large con¬ 
tribu tion from the gaseous Ca + in the ISM/C GM (see, 
e.g., iZhu fc Menardll2013aHMurga et al .1120151b 


and for emission, it is defined as 

fVax 

^emission = / [#( A ) _ f] dA , (4) 

J A m i n 

where the integration range (A m i n < A < A max ) encloses 
the absorption/emission profile. 

Throughout the paper, unless otherwise specified, we 
estimate the measurement uncertainties for a given sam¬ 
ple by bootstrapping (i.e., with replacement) 100 times. 

3.3. The Composite Continuum-normalized, Spectrum 

Figure [3] presents the median composite spectrum of 
the 8620 ELGs at 0.6 < z < 1.2 in the wavelength range 
2200 A < A < 2900 A. In Figure H we zoom in on the 
most prominent absorption and emission features for a 
more detailed illustration, and also include the part in 
the wavelength range 3000 A < A < 4000 A. We have 
omitted the part between 2900 A and 3000 A since it is 
featureless (see Figure [2]). To guide the eye, we mark the 
rest-frame positions of the identified features. 

The line features in the NUV can be categorized into 
three primary group^HI based upon their origins: stel¬ 
lar photospheric absorption lines, nebular emission lines, 
and absorption/emission lines due to gas in the ISM and 
CGM. The last category includes combined effects caused 
by both the ISM and CGM and is the focus of this work. 
We briefly discuss the observed features in these groups 
before investigating the ISM/CGM features in detail. 

• Stellar photospheric absorption features 


• Nebular emission features 

Between 2200 and 2900 A, we identify two nebular 
emission features with high confidence: semi-forbidden 
C II] at about 2326 A and forbidden [O II] at about 
2471 A. The C II] feature is a blend including five transi¬ 
tions from the doublet and the triplet between the ground 
state and first excited state: C II] AA2324,2325 and 
C II] AAA2326,2327,2329. The [O II] feature is a dou¬ 
blet at 2470.90 A and 2471.10 A, due to the transitions 
between the ground state and the second excited state of 
0 + . For comparison, the [O II] AA3727,3730 doublet is 
due to the transitions between the ground state and the 
first excited state of 0 + . Both the C II] and [O II] emis¬ 
sion lines were observed in H II regions by IUE (iDufour 
[1987 1. As the low-ionization nebular emission lines in 
the optical spectra of SFGs, e.g., [O II] AA3727,3730, 
[N II] AA6548,6583 and [S II] AA6718,6733, they must 
be dominated by the emission from H II regions, where 
the UV photons from O/B stars ionized the elements 
in the surround ing gas to low-ionization states (e.g., 
lStromgrenlll939fl . though diffuse ioni zed gas in the ISM 
e.g., t he warm ionized medium fWIM. lMcKee fc Ostriken 
19771 . als o cont ributes to the integrated emission (e.g., 
Reynolds! Il98l . C+ can also be abundant in a 
cooler medi um (e.g ., photodissociation regions, or PDRs, 
iTielens fc Hollenbachl 11985] ) because the ionization po¬ 
tential of neutral carbon (11.26 eV) is lower than that of 
hydrogen or oxygen (both ~ 13.6 eV, see Appendix [All. 
It is interesting to note that, though C II] in the NUV 
has been little studied, the fine-structure emission of the 
C + ground term at 157.7 //in in the IR is known to be 


17 In the FUV, stellar wind features comprise another major 
class of spectral features. 


18 http: //www. inaoep.mx/~modelos/uvblue/uvblue .html 
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Blue: Emission-line Galaxies (ELGs) at 0.6<z<1.2 



Fig. 3.— The composite continuum-normalized spectrum of 8620 ELGs at 0.6 < z < 1.2 (blue). We show the rest-frame positions of 
emission features with vertical green dashed lines and those of ISM/CGM absorption lines with black dashed lines. The brown solid line 
marks the stellar photospheric absorption feature C III A2298. The red line shows the composite spectrum of 2310 quasar intervening 
absorption-line systems with VF^ 2796 > 2 A in the same redshift range. 


a major coolant in the ISM (iDalearno fc McCravlll972ll . 
[C II] A157.7 /im can be observed fo r objects at 2 ~ 1 in 
the submillimeter (e.g., iStacev et alJl2010 .). thus for the 
ELG targets in eBOSS and future BAO surveys. 

We tentatively identify [Ne IV] (not labeled) at about 
2424 A, a doublet at 2422.56 A and 2425.14 A, due to 
the transitions between the ground state and the first 
excited state of the triply-ionized Ne 3+ . In SFGs, 
this doublet is observed mo stly in planetary nebulae 
('e.g.. [Koenpen fc Allertl 19871 ). supernova remnants (e.g., 
iBlair fc Panagial 119871) . and other environments hotter 
than H II regions as the ionization potential of Ne 2+ 
(63.4eV) is larger than He + (54.4eV). 

At 3000 < A < 4000 A, besides Balmer recombination 
lines and [O II] AA3727,3730, we also observe [Ne III] at 
A = 3869.86 A. Highly sensitive to the ionization param¬ 
eter, the [Ne III] emission can be combined with [O II] 
and used to probe the metallic i ty and other properties 
of the ISM dNagao et al.l 120061 iLevesaue fc Richardsonl 
12014 1. There is also a He I line at 3889.74 A (not labeled), 
which is blended with Balmer H£. Between 2900 A and 
the Balmer break at 3646 A, we only detect the weak 
He I emission at 3188.67 A, which likely sits on top of 
He I abs orpt ion com monly associated with O/B stars 
(e.g., lMorrisonlll975ll . 

• ISM/CGM absorption and emission features 

Most of the absorption lines between 2200 and 2900 A 
are induced by gas in the ISM and/or the CGM, pro¬ 
viding us with a unique tool to probe the diffuse gas 
otherwise elusive. From low to high energy, we identify 
the following lines due to different species: 

• Mg IA2853 (UV1); 

• the Mg IIAA2796,2804 doublet (UV1 J3 

19 Compared to the absorption induced by gas in the ISM/CGM, 
we expect the intrinsic Mg IA2853 and Mg IIAA2796, 2804 absorp- 


• the Mn IIAAA2577,2594,2606 triplet (UV1); 

• Fe IIA2600 and A2587 (UV1), Fe IIA2383 and 
A2374 (UV2), Fe IIA2344 (UV3), Fe IIA2261 
(UV4), and Fe IIA2250 (UV5). 

Between 3000 and 4000 A, the Ca II H & K lines also 
trace a significant fraction of gas in the ISM/CGM. All 
these absorption lines are also commonly seen in inter¬ 
vening quasar absorption-line systems (see next subsec¬ 
tion). 

We identify four non-resonant Fe II* emission lines 
with high confidence: 

• Fe II* A2626 and A2613 (UV1), Fe II* A2396 (UV2), 
and Fe II* A2366 (UV3). 

To illustrate the relationships between the resonant ab¬ 
sorption and these non-resonant emission lines, we use 
the Fe II UV1 multiplet group as an example. 

Figure [5] presents the energy-level diagram of Fe II 
UV1, showing the transitions between the ground state 
and the first excited state of Fe II. We refer the reader 
to Appendix [A] for a detailed description of the sym¬ 
bols and terms. Fe IIA2600 is the transition between 
the lowest (ground) energy level (with J = 9/2) of the 
ground state and the lowest level (with J = 9/2) of the 
excited state. Because the second lowest energy level of 
the ground state has the total angular momentum num¬ 
ber J = 7/2, the excited electron after the resonant ab¬ 
sorption has a high probability, i.e., a high Einstein A 
coefficient, to spontaneously decay to this level, releas¬ 
ing a fluorescent photon at a slightly longer wavelength 
A = 2626.45 A. The pair Fe IIA2586 and Fe II* A2613 
also belong to UV1, though with the second energy level 
of the excited state as the higher anchor level. 

The presence of these non-resonant emission lines im¬ 
ply that many of the resonant absorption lines above 
must be blended with emission filling in. For example, 

tion in the s pectra of O / B stars to be very weak, even for metal-rich 
stars (e.g., IRodnciiez-Menno et aflISOOa ). 
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A [A] 


Fig. 4.— The composite continuum-normalized spectrum of ELGs (blue), zoomed in on prominent features and compared with the 
composite spectrum of quasar absorbers (red) as in Figure [3] In addition, the bottom panel shows the spectra at 3000 A < A < 4000 A, 
where we have shifted the quasar absorber composite upward by 0.2 for clarity. 
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Fig. 5.— The energy diagram for the transitions in the first UV 
multiplet group (UV1) of Fe II, between the ground state and the 
first excited state. See Appendix A for a full description of the 
symbols and terms. 


the Einstein A coefficient for Fe IIA2600 (2.35 x 10 s s _1 ) 
is over six times higher than that for Fe II* A2626 
(3.53 x 10 7 s _1 ), so an excited electron has a higher 
probability to decay to the lowest level and release a 
Fe IIA2600 photon, though if the Fe II optical depth 
is high, the emitted photon will be absorbed again. The 
exact amount of emission infill depends on the optical 
depth of the relevant transitions and the balance between 
absorption and emission in the multiple-scattering pro¬ 
cess. 

We find the velocity profiles of both absorption and 
emission lines to be asymmetric, skewed towards nega¬ 
tive values, i.e., in the blueshift direction, indicating a 
larger fraction of the gas sources that cause these fea¬ 
tures are flowing outwards than inwards. The profiles 
of the emission lines appear to be similar, while the de¬ 
gree of asymmetry of the absorption profiles varies from 
line to line. For example, compared to Fe IIA2586 or 
Fe IIA2600 as shown in the second panel in Figure [4] 
the Mg IIAA2796, 2804 lines are bluesliifted from their 
rest-frame positions by a larger amount. Considering 
these species have similar ionization potentials (7.90 eV 
for neutral Fe and 7.64 eV for Mg) and likely co-exist 
spatially, the difference must be due to a larger amount 
of emission infill on top of the Mg II absorption. 

For a better understanding of the absorption and emis¬ 
sion in the composite spectrum of ELGs, we present a 
comparison with a composite spectrum of intervening 
quasar absorption-line systems in Section 13.41 and with 
a composite spectrum of local SF regions in Section 13.51 
We will then investigate the physical processes with a gas 
flow model in Section |4] 


The composite spectrum of ELGs exhibits the ab¬ 
sorption lines commonly seen in intervening quasar 
absorption-line systems. It is interesting to compare 
the so-called “down-the-barrel” spectra of ELGs with 
the intervening quasar absorber spectra. We select the 
absorbers from the JH U-SDSS metal absorber catalog 
dZhu fe Menardl 1201. 311). updat ed to the 12th Data Re¬ 
lease (DR12. lAlam et all l2015l'Pl We choose strong ab¬ 
sorbers with Mg IIA2796 rest equivalent width Wg 2796 > 
2 A, because there has been evidence showing that 
a large fraction of these strong absorbers are physi¬ 
cally asso ciated with the CGM of strong SF galax¬ 
ies (e.g.. iBergeron fo Boissel 1 19D lb iNorman et, al l 119961 : 
iBouclie et al.l 120071 : iNestor et all 1201 lb lLan et all 120111) . 
We select the 2310 such strong absorbers at 0.6 < z < 
1.2, the same redshift range as the ELGs and construct a 
median composite continuum-normalized spectrum. For 
more details regarding the quasar absorbers and how we 
estimate the continua of backgroun d quasars, we refer 
the reader to iZliu fc Menardl (|2013bH . 

We overplot the composite spectrum of quasar ab¬ 
sorbers in red in Figured and |4j Note the absorber spec¬ 
tra are based on luminous quasar spectra from SDSS I- 
III, so the S/N of their composite is orders-of-magnitude 
higher than the ELG spectra. When constructing the 
composite spectrum, we shifted the absorber spectra to 
the rest frame of the absorbers, so the absorption profiles 
are centered on their rest-frame positions. We find two 
major differences between the ELG composite and the 
quasar absorber composite: (1) There is no detectable 
non-resonant emission in the quasar absorber compos¬ 
ite, even though its S/N is orders-of-magnitude higher. 
(2) The ratios of the absorption lines are different. For 
example, the strength of Fe IIA2344 is similar in both 
composites, while Fe IIA2383 and Mg IIAA2796,2804 
are much stronger in the quasar absorber one. We note 
that the line ratios of quasar absorbers depend both on 
the strength (Wg 279e ) and redshift (e.g., Figure 3 in 
iDev et abl 1 201 all . but at Wg 2796 > 2 A and a given red¬ 
shift, the strength dependence is weak and has no effect 
on any of our conclusions. In addition to the two main 
differences, we do not detect Ti II absorption lines in the 
ELG composite, which we suspect is due to the limited 
sample size and the low S/N. 

In Section |4j we show that the line ratio difference is 
caused by the emission infill present in the ELG compos¬ 
ite but absent in the quasar absorber composite. 

3.5. Comparison with Local Star-forming Regions 

Spectroscopic observations of SFGs, or astronomi¬ 
cal sources in genera l, have been scarce in the NUV. 
iLeitherer et all (( 2011 4 compiled a UV spectroscopic at- 
la^P 7 ! of local SF galaxies observed with the Faint Ob¬ 
ject Spectrograph (FOS) and the Goddard High Res¬ 
olution Spectrograph (GHRS) on HST. The atlas in¬ 
cludes small-(physical) aperture spectra of 15 regions of 
nine SFGs with coverage between 2200 A and 3200 A, 
providing a rare opportunity for a direct compari¬ 
son of the integrated NUV SEDs at different phys¬ 
ical scales. From this compilation, we select nine 


3.4. Comparison with Intervening Absorbers 


20 http: //tow . pha. jhu. edu/~gz323/ j husdss 

21 http://www.stsci.edu/science/starburst/templ.html 
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Blue: Emission-line Galaxies (ELGs) at 0.6<z<1.2 (Aperture ~ 15 kpc) 



Fig. 6. — The composite continuum-normalize d spectrum of ELGs ( blue), the same as in Figure \3\ but compared with the composite 
spectrum of local star-forming regions (red) from [Leitherer et al.l (120111) . 


spectra of six galaxies with relatively high S/N and 
strong absorption: NGC1569, NGC2403 (all 3 spec¬ 
tra), NGC4569, NGC 5055, NGC 5253 (1 and 3), and 
NGC 5457 (NGC 5455). Note the different spectra of one 
galaxy are independent from each other, originating from 
different SF regions in that galaxy. Fo r deta ils re gard ing 
their atlas, we refer the reader to lLeitherer et aOl20111 1. 
With the selected nine spectra, we then construct the 
composite continuum-normalized spectrum following the 
same procedure as for the eBOSS ELGs. Before a careful 
comparison, we here emphasize several characteristics of 
the observations: 

• The aperture sizes of HST FOS/GHRS are one 
to a few arcseconds, and for all the nine spec¬ 
tra, correspond to physical sizes sm aller than 40 pc 
(2 — 37pc, see Table 6 in iLeitherer et al.1 120111) . 
They are nearly three orders-of-magnitude smaller 
compared to the aperture size for the eBOSS ELGs 
at 0.6 < 0 < 1.2, which is about 15 kpc (for 2 "0 

• The wavelength calibration of the FOS/GHRS 
spectra is largely based on the absorption lines in¬ 
duced by the gas in the ISM in the Milky Way, 
which are often blended with the lines of the low- 
redshift extragalactic sources. We find that, based 
on the position of C III A2298, the composite spec¬ 
trum is offset to longer wavelength by about 0.7 A 
(about gikms^ 1 ), we therefore correct the wave¬ 
length by —0.7 A. 

• The spectral resolutions of FOS/GHRS are one to 
a few hundred kms - , a factor of 2 — 4 lower than 
that of the BOSS spectrographs. 

The last two characteristics prevent us from making a 
quantitative comparison of the observed velocity profiles. 

In Figure [6l we compare the ELG composite (in blue) 
with the composite spectrum of the nine local SF re¬ 
gions (in red) at 2200 A < A < 2900 A. We present 

22 We note that the average seeing at the SDSS telescope is 
about 1.5”. 


a zoomed-in version in Appendix B. The two spectra 
share common absorption features across the wavelength 
range, including both the stellar photospheric absorp¬ 
tion line C III A2298 and the absorption lines due to the 
ISM/CGM. We also observe some common weak absorp¬ 
tion lines that must be intrinsic to the underlying stel¬ 
lar continuum. As discussed in Section 13.31 we do not 
identify these weak lines yet since they are still poorly- 
understood. We find the following main differences: (1) 
There is no detectable line emission, either the nebular 
lines (C II], [O II]) or the non-resonant lines (Fe II*), in 
the composite of local SF regions. (2) The absorption 
line ratios are different. For instance, the strength of 
Mg II AA2796, 2804 is similar in both composites, while 
Fe II A2586 and Fe II A2600 are about 50% weaker in the 
composite of local SF regions. We note that in the com¬ 
posite of local SF regions the ratios of the lines are sim¬ 
ilar to those in the quasar absorber composite, though 
the absolute absorption strength is about 50% weaker. 

The reasons for the non-detection of different emission 
lines are likely different. The [O II] doublet at 2471 A, 
mostly associated with H II regions, is very weak com¬ 
pared to its lower energy counterpart [O II] AA3727, 3730 
and may be buried in the noise. The C II] emission, 
which is strong in the ELG spectra, is more extended 
than H II regions (as observed through t he fine -structure 
emission [C II]A157.7^m, e.g.. iPineda et alJ 12013 1 be¬ 
cause of the lower ionization potential of neutral carbon, 
and the small (physical) apertures of FOS/GHRS did not 
capture enough C II] photons. The non-resonant Fe II* 
lines are not detected also because the emission is ex¬ 
tended, though the emission mechanism is different from 
C II]. 

As in the comparison with quasar absorbers, the dif¬ 
ference in absorption line ratios must be due to emission 
infill in the ELG spectra. We present a more detailed 
discussion below in the context of a gas outflow model. 

4. INTERPRETATION 

We have shown that the NUV composite spectrum of 
eBOSS ELGs display preferentially blueshifted absorp¬ 
tion, induced by neutral and singly-ionized species, Mg I, 
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Fig. 7.— The spherical outflow model. The color scale indicates the mean line-of-sight velocity of emission/absorption at a given 
position, assuming velocity v(r) oc r“, where a = 2 is arbitrarily chosen. Resonant absorption takes place in front of the background 
light source (stars). Fluorescent photons, resonant or non-resonant, are scattered isotropically and only those scattered into the line of 
sight can be captured. The aperture size of the eBOSS fibers is 2", corresponding to about 15 kpc at 0.6 < z < 1.2. The aperture 
sizes of 1 1ST F OS/ GHRS are one to a few arcseconds, corresponding to less than 40 pc for the local star-forming regions observed. See 
IScarlata K: Panagial 120151 1 for a similar model. 


Mg II, and Fe II. In addition, we detected non-resonant 
Fe II* emission lines, which also exhibit a preferentially 
bluesliifted profile and are not detected in either quasar 
absorber spectra or small-aperture spectra of local SF 
regions. These observed properties indicate that the gas 
causing the absorption and emission is predominantly 
flowing outwards, and the outflows must be driven by 
the strong SF activities in the ELGs and extend to large 
galactic scales. 

Galactic-scale outflows driven by star formation have 
been observed for over half a century, such as from 
the starburst galaxy M82 (e.g., iLvnds fe; Sandagell 19631 
I Bland fe Tull vl fl988h . The physic s of galactic winds has 
been extensively studied (e.g., iHeckman et al.l I1990L 
2000), though it is not yet conclusive due to the com¬ 
plexity of baryon processes involved. For our purposes, 
we circumvent some of the complex processes, such as 
the origins and properties of the wind and gas, instead 
we introduce a phenomenological gas outflow model and 
interpret our data with an observation-driven approach 
in the context of the model. 

4.1. A Spherical Gas Outflow Model 

We describe the model in three steps. First, we present 
the key geometrical and physical characteristics of the 
model. Second, we consider the properties of the obser¬ 
vations of the integrated spectra along the line of sight. 
Finally, we discuss quantitatively the radiative transfer 
processes and the model predictions for the SFG obser¬ 
vations. 

• Basics of the gas outflow model 

Because of the statistical fashion of our composite anal¬ 
ysis, we construct our model for the average of an ensem¬ 
ble of galaxies, not for a single source. We illustrate the 
model in Figure rfPl and first emphasize the following two 
key characteristics. 

23 The image in the center is a composite image of the central re¬ 
gion of M82 with different orientations. The original image is from 
http://hubblesite.org/gallery/album/galaxy/pr2006014a/, 


• Spherical symmetry A basic outflow model for an 
individual galaxy includes a density profile n(r) (in num¬ 
ber) and a velocity profile v(r ), both as a function of 
the vector position r. Observations have shown that 
the large-scale outflows driven by star formation are of¬ 
ten b ipolar, in the form of an expanding envelope (e.g., 
IHeckman et al .Ifl 9900 . In our composite analysis, we are 
averaging over all orientations randomly distributed and 
it is reasonable to assume spherically symmetric profiles 
n(r) and v(r). 

• Velocity distribution - Also due to the statistical na¬ 
ture of our approach, at a given galactocentric distance, 
there is a distribution of velocities. If gas accretion takes 
place around the ELGs and if the infalling gas is enriched 
with the species we are interested in, it will also induce 
absorption and contributes to the statistical signatures. 
At small scales, the gas in the ISM also contributes to 
the absorption, and its motions, ordered or disordered, 
also affect the observed velocity distribution. 

A statistical gas flow model therefore includes an aver¬ 
age density profile, n(r), an average velocity profile, v[r ), 
and a velocity dispersion profile, cr(r). We expect the di¬ 
rection of the average velocity v(r ) to be outwards as we 
expect there is more gas flowing outwards than inwards. 
The velocity dispersion then accounts for the contribu¬ 
tions at different velocities from outflows, inflows, and 
motions of ISM at small scales. In addition, in obser¬ 
vation, we need to consider the finite instrumental res¬ 
olution and redshift precision, which can be effectively 
included in the velocity dispersion through convolution. 
For the eBOSS ELGs, the mean spectral resolution is 
about {TV) ~ 2000 (60 - 70kms _1 ) and the mean red- 
shift precision is about 20kms^ 1 at redshift z ~ 0.8. 

Our model is in princ i ple si milar to the one introduced 
bv IScarlata fe Panagial (I2015L see also iRubin et al.ll201 ll : 
iProchaska et al.l 120111) . which the authors used as the 
basis in their radiative transfer simulations to interpret 
space-based observations in the FUV, though here we 

courtesy of NASA, ESA, and The Hubble Heritage Team 
(STScI/AURA). 
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emphasize the statistical aspect of our model. Consid¬ 
ering the sample size and S/N of our current data, we 
cannot yet place strong constraints on the model details, 
e.g, the functional form of the profiles or the parameters. 
In Figure 0 for illustration purposes, we assume a power 
law for the average velocity profile v(r) ~ r a with a = 2, 
which is arbitrarily chosen. We leave the full modeling 
to future work, and focus on the general properties of the 
model predictions below. 

• General model predictions 

Based on the model, we can predict the following gen¬ 
eral properties of the absorption and emission features in 
the integrated spectra along the line of sight. 

i. Origins and aperture dependence - In the model, 
the observed absorption and emission have different 
origins. The absorption is induced by gas in front 
of the background light source, e.g., the stellar pop¬ 
ulations inside the galaxy. The re-emitted (fluores¬ 
cent) photons, resonant or non-resonant, are scat¬ 
tered isotropically, so the observed emission origi¬ 
nates from gas located everywhere within the aper¬ 
ture, except from behind the galaxy due to occul- 
tation (see below). The strength of the emission 
therefore increases with the aperture size until the 
aperture encloses all the absorbing gas, while the 
strength of absorption depends little on the aper¬ 
ture size unless the column density varies signifi¬ 
cantly across the galaxy. 

ii. Net effect - The sum of absorption and emission in 
a given set of transitions, including all the resonant 
and non-resonant channels, would be zero if and 
only if the observer could collect all the re-emitted 
photons scattered into the line of sight with a very 
large aperture. However, because of the finite size 
of a galaxy, the photons behind the galaxy can¬ 
not penetrate the high -densi ty regi ons o f the galaxy 
to be captured. As lScarlata fe Panagial (12015T) . we 
refer to this effect as occultation. The finite size 
of galaxies makes the outflow model more compli¬ 
cated than that for stellar winds in which the star 
can be considered as a point source. The net ef¬ 
fect of absorption and emission is therefore always 
absorption. At redshift z ~ 0.8, the typical effec¬ 
tive radius of S FGs with M* ~ 10 1[ * is about 
2 — 4 kpc (e.g., [Williams e t all 1201(1 iWuvts et all 

usa. 

iii. Velocity profiles - If outflows are dominant (com¬ 
pared to inflows), the observed emission profile is 
asymmetric and preferentially blueshifted due to 
occultation of the redsliifted emission behind the 
galaxy. The profile of the absorption (even without 
emission infill) is also blueshifted, since it is only 
induced by the gas in front of the galaxy’s stellar 
populations (the light source). We expect the de¬ 
gree of blueshift is smaller for the emission profiles 
than absorption, because the re-emitted photons 
are scattered isotropically and only a negligible 
fraction of photons originating from the observed 
absorption are scattered into the line of sight, while 


other re-emitted photons within the aperture are 
less blueshifted. 

iv-1. Emission infill Like emission via the non-resonant 
channels, the model also predicts re-emitted pho¬ 
tons via the resonant channels, producing emission 
filling in on top of the absorption. The emission 
infill is not sufficient to compensate for all the ab¬ 
sorption so an observer always sees absorption (see 
point ii). However, because the emission profile is 
less blueshifted than absorption (point iii), the ob¬ 
served absorption profile is more blueshifted than 
the “true” absorption profile before emission infill. 
If the emission and absorption profiles are signif¬ 
icantly different, e.g., due to large outflow veloci¬ 
ties or large aperture (relative to the occultation), 
we also expect to see P-cygni-like profiles. The 
amount of emission infill depends on the aperture 
size, the galaxy (occultation) size, the permitted 
channels and their transition probabilities, and the 
optical depth, which determines whether the ob¬ 
served emission originates from a single or multiple 
scattering process and the relative strength of dif¬ 
ferent channels. The degree of emission infill (rel¬ 
ative to the absorption) is therefore different from 
line to line. Below we discuss this quantitatively 
for the lines in the NUV. 


• General radiative transfer processes 


In an expanding envelope, if the velocity gradient is 
large, the radiative transfer proce sses can be tre a ted un¬ 
der th e Sobolev app roximation dSobolevl fl96~dt iCastor 
119701: iRvbicki fc Hummed 11978 ). IScarlata fc Panagia 


( 20151 see also iProchaska et al. 1201 111 included an ex- 
cellent discussion of the general radiat ive transfer pro¬ 
cesses involved in a galactic outflow, and Prochask a et al.1 
(l20lllf also pioneered theoretical considerations of the 
NUV transitions. We refer the reader to those papers for 
more details. We here present a summary of the formulae 
most relevant to our analysis. 


• Absorption: In the Sobolev approximation, at a given 
position, the optical depth is given by 


r(r) 


— fiuM u ni(r) 

m e c 


du(r) 1 
dr 


(5) 


which is proportional to the density ni(r) at the lower 
level, the rest-frame wavelength of the transition A i u = 
X u i, oscillator strength f) u . and the inverse of the veloc¬ 
ity gradient (i.e., the thickness of the thin shell with the 
same velocity). We have ignored stimulated emission and 
angular dependence. Along the line of sight, the optical 
depth above applies to the velocity v(r), or equivalently, 
the wavelength A r = A; u [l — v(r)/c], at which the ab¬ 
sorption is given by R r = e~ T ^ r \ 


• Emission: In a single-scattering event, after the absorp¬ 
tion of a photon, the probability that the excited electron 
can decay from the upper level (it) to a given lower level 
(Z) is given by 


A u i 

Si A U i 


Pul = 


( 6 ) 
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where A U i is the spontaneous emission coefficient from 
the upper level u to the lower level i, and the summation 
is over all the possible channels in the lower state. The 
above equation ignores stimulated emission, collisional 
(de-)excitation, and also fine-structure emission within 
the same state. 

If the electron decays to the original level, in our case, 
the lowest level in the lower state, the re-emitted photon 
can be absorbed again, resulting in a multiple-scattering 
process. The escape probability of a resonant ph oton 
from a shell of optical depth t(v) is given by (|Mathisl 

mu 



and the fraction of the absorbed photons that are even¬ 
tually re-emitted via a non-resonant (nr) channel to a 
lower level l is 

OO 

fnr,l(T~) = Pnr,l ^ ( [7V(1 /?esc)] 

71=0 

_ _ Pnr,l _ 

l-p r (l- 0 „ c )’ W 

where p nT: i and p r are the probabilities of decaying to the 
non-resonant lower level l and the resonant lowest level, 
respectively (Eq. [6J , and we have omitted the upper level 
symbol u for simplicity. The fraction of the absorbed 
photons that are eventually re-emitted via the resonant 
channel and escape from the shell is 


/t ’( 7 ")—PrPesc ^ ( [Pr(l fiesc)\ 
n —0 

_ Prfiesc 

1 Pr(l esc) 


(9) 


When the optical depth of a given shell r( v) is shallow 
and the escape probability /3 esc ~ 1, we reach the single¬ 
scattering approximation, i.e., Eq. [6] 


fnr,l ~ Pnr,l ; and 

fr~Pr ■ (10) 

When the optical depth t(v) is deep and the escape 
probability /3 esc ss 0, all the re-emitted photons will 
be through the non-resonant channels, if there are any, 
with: 

r Pnr,l i 

fnr,i ~ ---and 

1 - Pr 

0 . ( 11 ) 


Note in both cases, summing over all channels gives 


E . fnr.i fr 1 • 
m.r.% 


( 12 ) 


• Radiative transfer processes in the NUV 

With the formalism above, we now investigate quan¬ 
titatively the absorption and emission lines in the NUV 
and their correlations in the context of the outflow model. 
We focus on eight resonant absorption lines, among 
which four have non-resonant channels: 

• Fe IIA2600 with Fe II* A2626 (UV1), 


• Fe IIA2587 with Fe II* A2613 and A2632 (UV1), 

• Fe IIA2374 with Fe II* A2396 (UV2), 


• Fe IIA2344 with Fe II* A2366 and A2381 (UV3), 
and the other four do not: 

• Mg IA2853, Mg IIA2804, Mg IIA2796, and 
Fe IIA2383 (UV2). 

We do not consider Fe IIA2261 (UV4), Fe IIA2250 
(UV5), and Mn IIAAA2577,2594,2606 here because of 
their lower S/N in the data. 

We present the relevant atomic data in Appendix [A] 
We note that we do not detect Fe II* A2632 from 
UV1, whose Einstein A coefficient is about half that of 
Fe II*A2613, and Fe II* A2381 from UV3 has Einstein A 
coefficient about half that of Fe II* A2366 and is blended 
with Fe IIA2383 from UV2. 

Within the context of the outflow model, we can now 
estimate the degree of the emission infill effect, i.e., the 
ratio of emission to absorption in the observed spectra, 
for the eight absorption lines in the NUV. 

iv-2. Degree of emission infill - The effect of emission 
infill depends on two main factors. For those with 
non-resonant channels, it depends on the fraction 
of resonant emission (/ r ). In the single-scattering 
approximation, we have f r ~ p r (Eqs. [TO] and 0 
and with the atomic data from Appendix S we 
have 


(13) 


When multiple-scattering events are considered 
(Eq. 0 , this order does not change though the rel¬ 
ative difference is smaller. We expect the effect of 
the emission infill to follow the same order. 


When there is no non-resonant channel, it mainly 
depends on the degree of saturation - the more 
saturated the absorption line is, the larger an ef¬ 
fect the emission infill has on the observation. Note 
the observed emission and absorption have different 
origins (point i). Based up on the elemental abun¬ 
dance (lAsplund et al.ll2009l) and oscillator strength 
(Appendix El , the degree of saturation should be 
in the order of absorption strength as 


tt/A2853 / tt/A2383 / tt/A2804 / tt/A2796 
vv Mgl < vv FeIl ^Mgll vv TV 


/ Mg II 


(14) 


Finally, we expect the degree of emission infill, 
manifested by the degree of blueshift (due to the 
difference in the profiles of emission and absorption 
(point iii)) and the change in the observed absorp¬ 
tion strength, to follow the same order given by 
Eos. fl3landfT4l 

For a direct comparison of observations with the model 
predictions, it is necessary to understand what the true 
absorption profiles (before the emission infill) are. In the 
next section tPD , we introduce an observation-driven 
method to reveal the true absorption profiles. 


4.2. Revealing the True Absorption P7vfiles 

To reveal the true absorption profiles, we make two 
assumptions: 
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Fig. 8.— The velocity profiles of the non-resonant emission lines, 
shown in the rest frame of the galaxies. Top panel : The observed 
profiles. Middle panel : The observed profiles normalized to the 
same amplitude. Bottom panel: The mean unified emission profile. 
The shaded region indicates the lcr uncertainties determined by 
bootstrapping. 

1. All the emission lines share the same normalized 
velocity profile. 

2. All the absorption lines share the same true nor¬ 
malized velocity profile, i.e., prior to the emission 
infill. 

In both cases, the lines are normalized to have the 
same amplitude. We make these assumptions accord¬ 
ing to our understanding of quasar absorption-line sys¬ 
tems from high S/N high-resolution spectroscopic ob¬ 
servations, which show that Fe II and Mg II (as well 
as Mg I when detect ed) usually trace each other (e.g., 
iChurchill et al .1 [20001) . We consider it reasonable to ex¬ 
trapolate these results to galaxy absorption lines. Since 
we work on the observed profile, i.e, i?(A) = e~ T ( x \ but 
not the optical depth r(A), saturation plays an impor¬ 
tant role and we will treat the Mg II lines separately in 
our method. 

Under the assumptions, our observation-driven 
method consists of two steps. 

1. We determine the common emission profile from 
the four non-resonant emission lines, which we call 
the unified profile and present in Section ld.2.11 

2. With the unified emission profile, we determine the 
unified true absorption profile with an iterative ap¬ 
proach. We describe this process in detail in Sec¬ 
tion [1221 

4.2.1. The Unified Emission Line Profile 

We combine the four non-resonant emission lines, 
Fe II* A2626, A2613, A2396, and A2366, to determine 



Fig. 9.- The asymmetry of the unified emission profile. Upper 
panel: The blue line and the shaded region show the mean unified 
emission profile and uncertainties as in the bottom panel of Fig¬ 
ure [H] The red dashed line shows a symmetric profile assuming the 
blueshifted side mirrors the observed redshifted side. Lower panel: 
The emission excess on the blueshifted side. The errorbar at the 
top right indicates the mean uncertainty at a given pixel (velocity). 

the unified profile, as shown in Figure [HI The top panel 
presents the observed velocity profiles, shown in the rest 
frame of the galaxies. In the middle panel, we nor¬ 
malize all the lines to have the same amplitude, and 
in the bottom panel, we present the mean normalized 
profile as the estimate of the unified emission profile. 
We also plot the lcr uncertainties determined from boot¬ 
strapping. We note that, because the far blue side of 
Fe II* A2613 at v < —600 km s^ 1 overlaps with the red 
side of Mn IIA2606, to avoid the contamination, we do 
not include the blue side of Fe II* A2613 at v < Okrns^ 1 
while calculating the mean profile. 

The unified emission profile appears to be asymmetric. 
We investigate its asymmetry further in Figure [9] In the 
upper panel, on top of the unified profile, we overlay a 
symmetric profile assuming the blue side mirrors the red 
side. We subtract the symmetric profile from the orig¬ 
inal one and show the result in the lower panel. The 
unified non-resonant Fe II* emission exhibits an excess 
on the blue side, with a confidence level higher than 2cr 
when integrated over —500 km s -1 < v < Okrns -1 , in¬ 
dicating that a larger fraction of emission is blueshifted 
than redshifted. 

4.2.2. The Unified Absorption Line Profile 

With the unified emission line profile, we determine the 
emission infill and the unified true absorption line profile 
simultaneously with an iterative approach. 

To proceed, we first investigate the observed velocity 
profiles in more detail. In Figure [TUI we present the ob¬ 
served profiles of the eight absorption lines in the rest 
frame of the galaxies, normalized to have the same ampli¬ 
tude integrated over —700 km s^ 1 <v< 300 km s . We 
have ordered the lines according to Eas. [T3land[Tdl with 
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Fig. 10.— The observed velocity profiles of the absorption lin es, prese nted in the rest frame of the galaxies and normalized to the same 
amplitude. The color scales indicate the order given by Eqs. I13l and ll4l 
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Fig. 11.— Left panels : Examples of emission-infill correction. The blue lines show the observed absorption velocity profiles, red the 
emission-corrected, and green the subtracted emission. The vertical dotted lines mark the rest-frame positions of the lines and the zero 
velocity corresponds to the wavelength of the one with lower energy in each panel. Right panels: Top - The emission-corrected absorption 
velocity profiles. Middle - The emission-corrected profiles normalized to the same amplitude. Bottom - The unified emission-corrected 
absorptio n p rofile. The shaded areas indicate the la bootstrapping uncertainties. The color scales in the top two panels are the same as 
in Figure FTOl 


bluer color indicating a larger predicted effect from emis- 
sion infill. Figure [TO] shows that: (1) the observed ab¬ 
sorption profiles are different; (2) the degree of blueshift 
follows the order predicted by the model, and (3) the 
most blueshifted line, Mg IIA2796, is about 200kms~ 1 
more blueshifted than the least blueshifted Fe IIA2374. 

Our iterative approach to determining the true absorp¬ 
tion profile consists of the following steps. 

1. We first use the normalized profile of Fe IIA2374 
as the initial guess of the unified true absorption 
profile, because the fluorescent emission after the 
Fe IIA2374 absorption is dominated by the non¬ 


resonant channel Fe II*A2396 lEo. Till . 

2. With the unified absorption profile estimated from 
the previous step, we fit for the amount of emis¬ 
sion that needs to be subtracted from the observed 
profile. More specifically, we express the observed 
profile (A) by 

R±(\) = R%£(\) + [R em i(A)-l] 

= {1 - a [1 - A)] } + b [i?em\(A) - 1] , 

(15) 

where £). e (A) and i? em i(A) are the unnormalized 
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Lines Lines 


Fig. 12.— The comparison of line ratios of ELGs with strong quasar absorbers (QSOABSs, left) and local SF regions (right). The blue 
circles show the obser ved line rat ios, and the green squares the emission-corrected ones. We use Fe IIA2374 as the anchor and order the 
lines according to Eq. 1 13 1 and 1141 The errorbars indicate uncertainties in ELG line ratio measurements determined by bootstrapping, not 
including the uncertainties in the measurements of QSOABSs or local SF regions. 


true absorption and emission profiles, respectively, 
and i?"£sW an d -ft™ 1 ;(A) are the unified normal¬ 
ized absorption profile from the previous step and 
the unified emission profile from Section 14.2.11 re¬ 
spectively. We perform a least-squares fit for the 
coefficients a and b. 

3. We normalize the new absorption profiles R^sW 
from the fitting in Step 2, calculate the mean as 
the new estimate of the unified absorption profile, 
and then repeat Step 2. 

As discussed above, saturation requires special attention 
when it is severe. We set R at saturated pixels to 
be zero and do not include Mg II while estimating the 
unified profile in Step 3. We iterate the steps until the 
unified absorption profile and the coefficients a and b 
converge. In practice, we find that three iterations are 
sufficient to reach convergence. 

We show the results in Figure [TlJ On the left, we 
show examples of the decomposition fEa. fTSl) . with the 
emission-infill indicated by the green dashed lines. The 
emission-corrected absorption profiles, shown with the 
red lines, are deeper but less blueshifted than the ob¬ 
served ones (blue). On the right, in the top panel, we 
show the emission-corrected profiles of all the eight ab¬ 
sorption lines, with color scales the same as in Figure [TUI 
The middle panel shows the emission-corrected profiles 
normalized to the same amplitude, which we use to esti¬ 
mate the unified true absorption profile. We present the 
unified profile in the bottom panel, together with the 
uncertainties estimated via bootstrapping. Note that af¬ 
ter emission-infill correction, the Mg II lines are heavily 
saturated at the fine centers and are not used in the cal¬ 
culation of the mean profile. 

The unified absorption profile is apparently asymmet¬ 
ric and preferentially blueshifted. With the emission- 
infill correction and the true absorption profiles esti¬ 
mated, we now investigate the details of the observed 
and emission-corrected profiles. 


4.2.3. Non-parametric characterization of the true 
absorption profiles 

To characterize the absorption profiles, we choose to 
use non-parametric variables, line ratios based on rest 
equivalent width (Eq. [3]) and velocity offsets. 

• Line ratios 

Because of the different degrees of emission infill, the 
changes in the absorption strength vary from line to line. 
We compare line ratios before and after the emission-infill 
correction to study if the effect of emission infill follows 
the order of Eqs. [12] and HU 

We measure the rest equivalent width of the emission- 
corrected profiles by integrating over —700kms~ 1 < 
v < 300kms _1 . For the observed profiles of fines 
other than the Mg II doublet, we integrate over the 
same velocity range. We select the integration veloc¬ 
ity range —700kms _1 < v < 200 km s -1 for the ob¬ 
served Mg IIA2796 and — 600 km s -1 < v < 300 km s^ 1 
for Mg IIA2804 to avoid the contaminations from each 
other (see Figure flOl). 

To calculate fine ratios, we select Fe IIA2374 as the an¬ 
chor, i.e., the common denominator, which has negligible 
emission infill. In Figure 1121 we compare the line ratios 
before and after the emission-infill correction with the 
line ratios in the composite spectrum of strong quasar 
absorbers (QSOABSs, in the same redshift range) on the 
left and with those of local SF regions on the right. The 
variables plotted are the ratios of line ratios: 


/ W Xi \ 

\^w A2374 ) 


ELG 


/ W Xi \ 

1 ^/A2374 J 


QSOABS 


and 


W Xi A 
W x2374 / 


ELG 


W Xi 
W X2374 


Local SF 


( 16 ) 


where A,; represents the fine names. In the figure, we have 
ordered the lines according to the predictions of Eqs. [13] 
and[H 
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Lines 


Fig. 13.— The non-parametric characterization of the absorption 
velocity profiles. We show i^o {blue), V 50 {green) and V 20 (red), 
before {circles) and after {horizontal bands) the emission- infi ll cor¬ 
rection. The order of the lines is the same as in Figure Il2| The 
width of the horizontal bands indicate the uncertainties. All the 
uncertainties are determined by bootstrapping. 

Figure [12] shows that the observed line ratios in the 
composite spectrum of ELGs differ from those of strong 
quasar absorbers and also local SF regions, and the de¬ 
gree of the difference basically follows the predicted or¬ 
der. We find that the effect of emission infill (on line 
strength and ratio) can be larger than a factor of two, 
e.g., for Fe IIA2383 and Mg IIAA2796, 2804. After the 
emission-infill correction, we find that the line ratios are 
consistent in the spectra of all the sources. 

We stress that our method of estimating the emis¬ 
sion infill is observation-driven, completely independent 
of any model. We only made the two assumptions about 
the unified emission and absorption profiles as elaborated 
at the beginning of this section. We also did not use any 
information about the strong quasar absorbers or the lo¬ 
cal SF regions. The agreement of the final line ratios 
is therefore not by construction, but rather the result 
of the method. As we expect that the gases inducing 
the absorption lines in different sources have similar ori¬ 
gins, either supernova yields, stellar mass losses or other 
sources, and the line ratios in different source spectra 
should be in agreement, we consider our iterative fitting 
method successful in estimating the emission infill and 
determining the true absorption profiles. 

• Velocity offsets 

To characterize the velocity profile, we define a velocity 
offset variable u xx to be the velocity where a fraction of 
xx per cent of the absorption is at velocity v > v x m 
In Figure E] we show vgo, V 50 and V 20 of the observed 
absorption profiles and of the unified profile. The order 
of lines is the same as in Figure El given by Eos. Eland 

24 We note that this is a characterization of the total absorption 
profile, including contributions from both outflows and ISM. 



Fig. 14.— Upper panel : The comparison of the unified, emission- 
corrected absorption velocity profile {red), the unified emission pro¬ 
file {blue), and a Gaussian profile with the width of 108kms _1 
{green). For display purposes, we have flipped the absorption pro¬ 
file and also adjusted the profiles so that they have roughly the 
same peak value. The shaded regions show the lu bootstrapping 
uncertainties. Lower panel : The comparison of the normalized 
[O III] A5008 emission profile for ELGs (at 0.6 < z < 1.2, magenta) 
with the same Gaussian profile {green) as in the upper panel. The 
uncertainties of the [O III] profile are about the same size as the 
line width. 

El 

Figure [FJ] (see also Figure [TUI) shows that, for the ob¬ 
served profiles, the velocity offsets of different lines are 
different, and the order follows the one predicted by 
the radiative transfer considerations in Eqs. E] and El 
and the difference in the emission and absorption pro¬ 
files (point iii in Section 14.11) . The unified emission- 
corrected profile is still asymmetric, but with a smaller 
degree of blueshift: the 50% velocity offset H 50 is about 
—lOkms -1 , and |i> 8 o| is about 160kms -1 , larger than 
1 ^ 20 1 ~ llOkms , Our fitting result shows that with¬ 
out correcting for the emission infill, ugo, the 80% veloc¬ 
ity offset, can be overestimated by over a factor of two 
for the lines most severely affected, e.g., Fe IIA2383 and 
Mg IIAA2796,2804. 

4.3. Discussion 

We now compare the observations, including the re¬ 
sults of our emission-infill correction method, with the 
predictions of the spherical outflow model presented in 
Section l4~Tl We go over the predictions point by point. 

[i.] Aperture dependence - The physical aperture size 
of the eBOSS ELG spectra is about 15 kpc, while that 
of the HST FOS/GHRS spectra of the local SF regions 
is smaller than 40 pc. In the spectra of the local SF 
regions, we do not detect the non-resonant emission that 
otherwise persist in the ELG ones. This agrees with the 
model and is because the HST FOS/GHRS aperture size 
is too small to capture the extended fluorescent emission. 

[ii.] Net effect - The net effect is always absorption due 
to occultation in the model, even if the aperture encloses 
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all the emission scattered into the line of sight. In the 
data, when we sum up all the absorption and emission 
in a given set of channels, the net result is absorption. 


[iii.] Velocity profiles - The outflow model predicts that 
both the emission and (emission-infill corrected) absorp¬ 
tion profiles are blueshifted, with the absorption more 
so than the emission. Figures [9] and [13] quantify the 
asymmetry and the blueshift of the profiles individually, 
though in different ways. To compare the two profiles 
directly, we show them together in Figure 1141 For dis¬ 
play purposes, we have flipped the absorption profile and 
also normalized them so that they roughly have the same 
peak value. The absorption profile is more blueshifted 
than the emission, as in the model. 

As discussed in the basics of the model, the composite 
spectra include effects from not only the outflows, but 
also the inflows, the motions of the ISM in the galaxy, 
the instrumental resolution as well as the redshift pre¬ 
cision. If the aperture is large and collects all the re¬ 
emitted photons along the line of sight, and if the out¬ 
flowing gas is extended to much larger scale than the 
galaxy, we expect the red side of the emission profiles to 
be broader than nebular lines, extended to further red. 
On the other hand, if a substantial fraction of the gas 
is falling in onto the galaxy at high velocities, we also 
expect the red side of both the emission and absorption 
profiles to be broader. In the lower panel of Figure [T4l 
we show the [O III] A5008 profile of the NUV sample. We 
note that ELGs at z > 1 do not have [O III] coverage, 
but they account for a small fraction (10%) of the sam¬ 
ple. We find the [O III] A5008 profile is well-represented 
by a Gaussian profile and the best-fit width (a) is about 
108kms _1 . We overplot this Gaussian profile in both 
panels for comparison. 

The mean spectral resolution of the BOSS spectro¬ 
graphs is about 60 — 70kms~ 1 and the average redshift 
precision of the eBOSS ELGs at redshift 0.6 < z < 1.2 is 
about 20kms _ . The width of the nebular emission line 
profile (108 kms^ 1 ) must therefore be dominated by the 
intrinsic rotation and disordered motion of the ISM along 
the line of sight, which account for about SSkms^ 1 when 
we subtract the spectral resolution and redshift precision 
by quadrature. This is consistent with measurements 
of kinematic properties and the Tully- Fisher relation for 
bright galaxies at these redshifts (e.g, ., lYogt et aI1ll996t 
IWeiner et al.ll2006[ I Miller et alJl201lh . 


We find that, on the red side, both the emission and ab¬ 
sorption profiles are consistent with [O III] A5008 within 
the uncertainties. This suggests that, within the uncer¬ 
tainties of our data, we have not observed evidence for 
outflows on scales larger than the galaxy nor the evidence 
for inflows. The larger-scale outflows would extend the 
emission profile to higher velocities (farther on the red 
side), and the inflows would extend both the emission 
and absorption profile. However, the former (that we did 
not see evidence for larger-scale outflows) could be be¬ 
cause the aperture size (15 kpc) is not sufficiently large, 
while the latter could be becaus e the inflow velocities 
(e.g., < 200 km s -1 . iRubin et all 120121 ) are comparable 
to the ISM motions and the emission/absorption effects 
of the inflows and the ISM are blended together. 

If we assume the emission-corrected absorption on the 
red side solely originates from the ISM, we can decom¬ 


pose the total absorption into a part due to the ISM 
with a sym metric profile and the other due to the out¬ 
flows (e.g., (Weiner et al.ll2009rr4 . Applying this decom¬ 
position method to the unified absorption-corrected pro¬ 
file, we obtain a blueshifted excess profile, more extended 
than the emission excess at v < 0 km s _1 as shown in Fig¬ 
ure [H with maximum velocity about —600 km s -1 . The 
ratio of the excess to the subtracted symmetric profile, 
which represents the ratio of the amount of outflowing 
gas to that of the ISM assuming the decomposition is 
ideal, is about 1:3. However, because there could be a 
broad distribution of outflow velocities, i.e., a large cr(r) 
in the outflow model even without the ISM contribution, 
it is likely that the symmetric component also includes a 
large contribution from the outflowing and/or inflowing 
gas, in which case the contribution from the ISM is much 
smaller than 3/4. 

[iv.] (Degree of) Emission infill The model predicts 
that, in large-aperture spectra, there is emission filling 
in on top of the resonant absorption, and because the 
emission profile is less blueshifted, the infill results in an 
observed absorption profile that is more blueshifted. The 
degree of the emission infill depends on the transition 
probabilities of the permitted channels and the degree of 
saturation for those without non-resonant channels. We 
have demonstrated that the data agree with these pre¬ 
dictions in Figure ITOlfTSl In particular, we show that the 
blueshifts, line ratios and velocity offsets of the observed 
profiles follow the same order as predicted (Eqs. [13] and 
El). The observed uso, the 80% velocity offset, can be 
overestimated by a factor of two compared to that in the 
emission-corrected profile. After the emission correction, 
the line ratios in the spectra of ELGs are consistent with 
those of strong quasar absorbers and local SF regions. 

The different degre es of blueshift for diffe rent lines 
were also suggested bvlProchask a et al.l (120111 1 and were 
observ ed in the Keck spect ra of SFG s by Erb et al.1 
(12012L see their F igure 7). lErb et a.l~l 1)2012 . see also 
iKornei etTa l. 2013) also found a variety of Mg II profiles 
in their individual spectra, with some showing emission 
that might originate from H II regions. In our model, 
we have ignored such contribution. With our data, we 
cannot yet quantify the effect of the emission from H II 
regions on the line profiles in the composite analysis. 

We do not observe P-cygni-like profiles in the compos¬ 
ite spectrum of the full sample. This is likely because 
the difference of the emission and absorption profiles on 
the red side is not large (Figure [14]) and the amount of 
emission infill is not sufficient, since P-cygni-like profiles 
require a large amount of emission infill that is more 
extended on the red side. In the next section, when 
studying the [O II] AA3727,3730 dependence, we show 
that the Mg II absorption features P-cygni-like profiles 
for the subsample with the higher [O II] AA3727, 3730 rest 
equivalent width. 

In summary, we conclude that our statistical, spherical 
outflow model can simultaneously explain the multiple 
observed properties of emission and absorption features 
in the NUV. 


25 We note IWeiner et all d200fj i applied this method to the ob¬ 
served Mg II absorption as they did not consider the emission infill. 
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Fig. 15.— The dependences of the rest equivalent width of the non-resonant emission lines on the [O II] AA3727, 3730 rest equivalent 
width (left) and luminosity (riqht). Note both the x axes are in logarithmic scale. 
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Fig. 16. — The dependences of the observed absorption velocity 
profiles of the Mg II lines on the [O II] AA3727, 3730 rest equivalent 
width (upper panel ) and luminosity (lower panel). 

5. CORRELATIONS WITH [O II] AA3727, 3730 

Observations have shown that outflow properties, 
such as the velocity, depend on galaxy pr operties 
(e.g., iRunke et al.1 12005L iTremonti et al.ll2007li . From 
the eBOSS pilot observations, we can measure the 
[O II] AA3727,3730 properties of the ELGs. We here 
study the dependences of the emission and absorption 
lines in the NUV on the total rest equivalent width 
(^[onf) and luminosity (L$™) of the [O II] doublet. 

For each variable, we divide the NUV sample into two 
subsamples, split at the median values (^^Fprq 8 ) = 

51.4 A and ^log 10 Ajonf) = 41-6 dex). We then perform 

the same analysis as for the full sample, including con¬ 
structing the composite spectra, calculating the unified 
emission profiles, estimating the emission-infill and de¬ 


termining the true absorption profiles. We present some 
of the details in Appendix [Cl includi ng th e distributions 
of W'^jn] 8 and logio^nf (Figure ICT]). the observed 
emission/absorption profiles (Figures IC2I and [C3ll . and 
the (unnor maliz ed) emission-corrected absorption pro¬ 
files (Figure [CD). We here discuss in detail the emission 
strength, the observed absorption profiles of Mg II, the 
emission-corrected absorption strength, and the unified 
velocity profiles. 

• Emission strength - Figure [15] shows the depen¬ 
dences of the non-resonant emission strength on the 
[O II] AA3727,3730 rest equivalent width and luminos¬ 
ity. Note for comparison, we have added the data points 
measured from the full sample in solid symbols, which 
are correlated with the measurements based on the two 
subsamples. Within the range probed, we find that the 
emission strength (in rest equivalent width) scales almost 
linearly with IVj( 3 r 7 p 8 . In the right panel, we show that 
the emission equivalent width also positively depends on 
the [O II] luminosity, but to a lesser degree. 

• Observed absorption profiles - In Figure flfil we present 
the observed Mg II profiles, which show the strongest 
dep ende nce on [O II] among the absorption lines (Fig¬ 
ure IC3D . The correlation appears to be stronger with 
the rest equivalent width than with the luminosity. This 
is due to the stronger dependence of the emission in¬ 
fill on the rest equivalent width, as suggested by Fig¬ 
ure HU The observed Mg II profile for the subsample 
with higher Wjp has a P-cygni-like shape, indicating 
a large amount of emission infill. 

• Emission-corrected absorption strength Figure [TT] 

shows the dependences of the emission-corrected absorp¬ 
tion rest equivalent width on the [O II] properties. Ex¬ 
cept for the saturated Mg II lines, other lines are posi¬ 
tively correlated with both and n®, with the 
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[OII1 EW Dependence of Absorption 
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Fig. 17.— The dependences of the rest equivalent width of the emission-correc ted absorption lines on the [O II] AA372 7, 37 30 rest equivalent 
width (left) and luminosity (right). The color scales are the same as in Figure flUl based on the orders given by Eqs. Il3l and ll4l 
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Fig. 18. — The dependences of the unified emission profile (left) and the unified emission-corrected absorption profile (right) on the 
[O II] AA3727, 3730 rest equivalent width (upper panels) and luminosity (lower panels). 
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dependence tentatively stronger for the former. 

• Unified velocity profiles - Figure [18] presents the unified 
emission and absorption profiles as a function of [O II] 
rest equivalent width and luminosity. Within the uncer¬ 
tainties, for both profiles, we do not find a dependence 
on either W^Jf 8 or Ljonp although the P-cygni-like 
shape of the observed Mg II absorption for the subsam¬ 
ple with higher requires the emission profile to 

be more extended on the red side than the absorption, 
unlike for the main sample (Figure [TJ)- Larger samples 
in the future will help pin down these dependences with 
high S/N. 

Among all the correlations, we find that the strongest 
one is between the rest equivalent widths of non-resonant 
emission and [O II], which also results in the strong de¬ 
pendence of the observed absorption profiles, especially 
of the Mg II lines, on the [O II] rest equivalent width. 
To the first order, the [O II] luminosity is an indicator of 
SFR, while the [O II] rest equivalent width is an indicator 
of specific SFR. Our results suggest that the properties 
of the emission and, to a lesser degree, the absorption are 
stronger functions of specific SFR than of SFR. However, 
considering the uncertainties due to our sample size, the 
exact correlations between the properties of the spectral 
features in the NUV and those of galaxies and their im¬ 
plications for galaxy evolution remain to be determined. 

6. SUMMARY 

The pilot observations of the emission-line galaxy 
(ELG) program in the Extended Baryon Oscillation 
Spectroscopic Survey (eBOSS) in SDSS-IV have obtained 
a sample of 8620 ELGs at 0.6 < 2 < 1.2, provid¬ 
ing a good opportunity for investigations of the near¬ 
ultraviolet (NUV) part of the spectral energy distribu¬ 
tions (SEDs) of star-forming galaxies (SFGs). We con¬ 
structed median composite continuum-normalized spec¬ 
tra to study the emission and absorption features in the 
NUV. Our main results are: 

• The median composite spectra of the ELGs feature 
non-resonant Fe II* emission and resonant absorp¬ 
tion due to Mg I, Mg II and Fe II. Both the emis¬ 
sion and absorption profiles are asymmetric, prefer¬ 
entially blueshifted, indicating ubiquitous outflows 
driven by star formation at 0.6 < z < 1.2. 

• We found a variety of velocity profiles for the ob¬ 
served absorption lines with different degrees of 
blueshift. 

• Comparing the ELG spectra with those of inter¬ 
vening quasar absorption-line systems in the same 
redshift range, we found they feature the same ab¬ 
sorption lines but with different line ratios. 

• We compared the eBOSS ELG spectra with the 
NUV spectra of the local star-forming regions taken 
with the Faint Object Spectrograph (FOS) and 
Goddard High-resolution Spectrograph (GHRS) on 
HST. The physical aperture size of the eBOSS ELG 
spectra at 0.6 < 2 < 1.2 is about 15 kpc, while the 
aperture size of the FOS/GHRS spectra of the lo¬ 
cal SF regions is less than 40 pc. We found the 
FOS/GHRS spectra also display the same (though 


weaker) absorption lines, but do not exhibit the 
non-resonant Fe II* emission. We also found differ¬ 
ent ratios for the resonant absorption lines. 

• We introduced a statistical, spherical outflow 
model, in which the observed non-resonant emis¬ 
sion is the fluorescent (re-emitted) photons after 
the occurrence of absorption that are scattered into 
the line of sight. The model predicts that there 
is scattered resonant emission Filing in on top of 
absorption, and the amount of emission infill de¬ 
pends on the transition probabilities of the allowed 
channels, resulting in the variety of the observed 
absorption profiles. 

• We developed an observation-driven, model- 
independent method to estimate the emission infill 
and reveal the true absorption profile. We showed 
that after the emission correction, the absorption 
line ratios in the ELG spectra are consistent with 
those in the spectra of strong quasar absorbers and 
local star-forming regions. 

• We demonstrated that the outflow model can ex¬ 
plain simultaneously the multiple observed proper¬ 
ties of the emission and absorption features in the 
NUV, including i) the aperture dependence, ii) the 
net effect, iii) the emission velocity profiles and the 
emission-infill corrected absorption profiles, and iv) 
the variety of the observed absorption profiles and 
the degree of emission infill. 

• Finally, we investigated the dependence of NUV 
features on the [O II] AA3727, 3730 rest equivalent 
width and luminosity and found that the strongest 
correlation is between the non-resonant emission 
strength (in rest equivalent width) and the [O II] 
rest equivalent width. 

Our observations provided strong evidence for ubiq¬ 
uitous galactic-scale outflows driven by star formation. 
Our analysis also demonstrated that the NUV window 
is an informative region in the spectrum. The se¬ 
ries of emission and absorption lines provides a new 
means to prob e the gas physics. The model we intro¬ 
duced (see alsolRubin et aLlfeoiUlProchaska et al.ll201 lb 
LScarlata fc Panagiall2015[T has many important implica¬ 

tions, such as the dependences of non-resonant emission 
on the aperture, outflow, and galaxy (occultation) sizes, 
and points the future investigations of outflow physics 
into new directions. For instance, it is of great inter¬ 
est to explore the surface brightness profile of the non¬ 
resonant emission, e.g., through narro w-band imag ing or 
spatially-resolved s pectroscopy (e.g., iRubin et al.ll201lL 
iMartin et al.1 [20131. to further study the scale depen¬ 
dence of outflows. 

Besides using the “down-the-barrel” spectra to probe 
the gas physics directly associated with the galaxies, we 
can also employ th e cross-correlation techniques deve l- 
oped recently (e.g., iSteidel et al.ll201(ilZhu et al .1120141 ). 
using absorption information induced in background 
source spectra, to probe the circumgalactic medium of 
foreground sources. Combining the two different types 
of observation will produce a more complete picture of 
the baryon processes in galaxy formation and evolution. 
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The sample size of the NUV spectroscopic datasets 
will grow by orders-of-magnitude in the next decade. At 
the conclusion of the ELG program, eBOSS will obtain 
spectra for about 200, 000 ELGs at z > 0.6, a sample 
over 2 0 times larger than t he one used in this paper. 
DESI dSchlegel et' al.l 1201 It iLevi et all l2013li and PFS 
dTakada et al.ll2014f) will obtain higher-resolution spec¬ 
tra with larger telescopes for about 20 million ELGs at 
higher redshift (z > 1), where more lines are redshifted 
into the optical. Based upon the details revealed in the 
composite spectra of less than 10, 000 galaxies with the 
2.5-meter SDSS telescope, we expect that NUV spec¬ 
troscopy will play an important role in future investiga¬ 
tions of the properties and evolution of galaxies. 
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TABLE Al 

The Ionization Potential" 


Atom 

N p 

An 

I-s-II 

II—>-111 

III-!-IV 

IV->V 

V->VI 

VI-)-VII 

VII-)-VIII 

H 

1 

0 

13.59843 







He 

2 

2 

24.58739 

54.41776 






Li 

3 

4 

5.391715 

75.64009 

122.4544 





Be 

4 

5 

9.322699 

18.21115 

153.8962 

217.7186 




B 

5 

6 

8.298019 

25.15483 

37.93058 

259.3715 

340.2260 



C 

6 

6 

11.26030 

24.3845 

47.88778 

64.49358 

392.0905 

489.9932 


N 

7 

7 

14.53413 

29.60125 

47.4453 

77.4735 

97.89013 

552.0673 

667.0461 

O 

8 

8 

13.61805 

35.12111 

54.93554 

77.4135 

113.8989 

138.1189 

739.3268 

F 

9 

10 

17.42282 

34.97081 

62.7080 

87.175 

114.249 

157.1631 

185.1868 

Ne 

10 

10 

21.56454 

40.96296 

63.4233 

97.1900 

126.247 

157.934 

207.271 

Na 

11 

12 

5.139077 

47.28636 

71.6200 

98.936 

138.404 

172.23 

208.504 

Mg 

12 

12 

7.646235 

15.03527 

80.1436 

109.2654 

141.33 

186.76 

225.02 

Al 

13 

14 

5.985768 

18.82855 

28.44764 

119.9924 

153.825 

190.49 

241.76 

Si 

14 

14 

8.151683 

16.34585 

33.493 

45.14179 

166.767 

205.267 

246.32 

P 

15 

16 

10.48669 

19.76949 

30.20264 

51.44387 

65.02511 

220.430 

263.57 

s 

16 

16 

10.36001 

23.33788 

34.86 

47.222 

72.5945 

88.0529 

280.954 

Cl 

17 

18 

12.96763 

23.81364 

39.80 

53.24 

67.68 

96.94 

114.2013 

Ar 

18 

22 

15.75961 

27.62967 

40.735 

59.58 

74.84 

91.29 

124.41 

K 

19 

20 

4.340664 

31.625 

45.8031 

60.917 

82.66 

99.44 

117.56 

Ca 

20 

20 

6.113155 

11.87172 

50.91315 

67.2732 

84.34 

108.78 

127.21 

Sc 

21 

24 

6.56149 

12.79977 

24.75684 

73.4894 

91.95 

110.68 

137.99 

Ti 

22 

26 

6.82812 

13.5755 

27.49171 

43.26717 

99.299 

119.533 

140.68 

V 

23 

28 

6.746187 

14.618 

29.3111 

46.709 

65.28165 

128.125 

150.72 

Cr 

24 

28 

6.76651 

16.48630 

30.959 

49.16 

69.46 

90.6349 

160.29 

Mn 

25 

30 

7.434018 

15.63999 

33.668 

51.2 

72.41 

95.604 

119.203 

Fe 

26 

30 

7.902468 

16.1992 

30.651 

54.91 

75.0 

98.985 

124.98 

Co 

27 

32 

7.88101 

17.0844 

33.50 

51.27 

79.50 

102.0 

128.9 

Ni 

28 

31 

7.639877 

18.16884 

35.187 

54.92 

76.06 

108.0 

132.0 

Cu 

29 

35 

7.72638 

20.29239 

36.841 

57.38 

79.8 

103.0 

139.0 

Zn 

30 

35 

9.394199 

17.96439 

39.7233 

59.573 

82.6 

108.0 

133.9 


" In unit of eV. Data are taken from the NIST-ASD database (IKramida et a l~] 1 2IJI lii . Only most abundant 
isotopes are included. 
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APPENDIX 
A. ATOMIC DATA 

The spectroscopic analysis presented in this paper requires accurate atomic data. In this appendix, we present the 
data we gathered for completeness and future references. 

We collect most of th e da ta from the National Institute of Standards and Technology Atomic Spectra Database 
(NIST-ASD. IKramida et aTll2014 ) and have cross-checked the information with either the original or new references 
we found in the literature. 

Most of the conventions used in spectroscopy are well established and we follow them whenever possible. However, 
we also found some ambiguous terms that m ay cause con f usion , mostly due to histo rical reasons . We try to be as 
specific as possible on those occasions. We use lMoore] (1 1 95(1 IT972I) . iCohen et al.1 (119871) . and lDrainel (120111) as the main 
references. We briefly summarize the nomenclature we adopt below. We only include those definitions that are most 
relevant to our discussion or may cause confusion. For a complete description, we refer the reader to the references 
above. We note our intention is not to give a comprehensive review or to resolve all the disagreements between the 
nomenclature in observational astronomy and that in quantum theory, but rather to make our presentation clear while 
maintaining consistency with the literature in astronomy. 

• Nomenclature 

o Roman numerals: To specify the spectrum of a z-fold ionized atom, we use a small capital roman numeral 
corresponding to z + 1 written with a thin space following the chemical symbol. For example, H I denotes the 
spectrum of neutral hydrogen (H°), Mg II that of singly-ionized magnesium (Mg + ), C III that of doubly-ionized 
carbon (C 2+ ). Recently, this convention has also been used by many authors to specify the ionized atom itself, 
not only the spectrum. 

o Electron configuration: An atomic electron configuration is indicated symbolically by 

(■ nl) k (n'l'f , (Al) 

in which k, k' , ... are the numbers of electrons with principal quantum numbers n, n 1 , ... and orbital angular 
momentum quantum numbers 1, l', ..., respectively. The orbital angular momentum numbers l = 0, 1, 2, 3, 
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TABLE A2 

The list of lines between 2200 A and 7500 A associated with star-forming galaxies 



Name 

Wavelength 
A (Vac.) 

A-ul 

s- 1 

Jlu 

Ei 

cm 1 

E u 

cm -1 

Multiplet 

No. 

ReP 

Fe IIA2250 

2249.88 

3.00E+06 

1.82E-03 

0 

44446.905 

UV5 

1,2 

Fe IIA2261 

2260.78 

3.18E+06 

2.44E-03 

0 

44232.540 

UV4 

1,2 

Fe II 

* A2270 

2269.52 

4.00E+05 

3.10E-04 

384.7872 

44446.905 

UV5 

1,2 

Fe II 

* A2281 

2280.62 

4.49E+06 

4.38E-03 

384.7872 

44232.540 

UV4 

1,2 

C III A2298 

2297.58 

1.38E+08 

1.82E-01 

102352.04 

145876.13 

UV8 

3 

C II 

A2324 

2324.21 



0 

43025.3 

UV0.01 

3,4 

C II 

A2325 

2325.40 



0 

43003.3 

UV0.01 

3,4 

C II 

A2326 

2326.11 



63.42 

43053.6 

UV0.01 

3,4 

C II 

A2328 

2327.64 



63.42 

43025.3 

UV0.01 

3,4 

C II 

A2329 

2328.83 



63.42 

43003.3 

UV0.01 

3,4 

Fe IIA2344 

2344.21 

1.73E+08 

1.14E-01 

0 

42658.244 

UV3 

1,2 

Fe II 

* A2366 

2365.55 

5.90E+07 

4.95E-02 

384.7872 

42658.244 

UV3 

1,2 

Fe IIA2374 

2374.46 

4.25E+07 

3.59E-02 

0 

42114.838 

UV2 

1,2 

Fe II 

* A2381 

2381.49 

3.10E+07 

3.51E-02 

667.6829 

42658.244 

UV3 

1,2 

Fe IIA2383 

2382.76 

3.13E+08 

3.20E-01 

0 

41968.070 

UV2 

1,2 

Fe II 

* A2396 

2396.35 

2.59E+08 

2.79E-01 

384.7872 

42114.838 

UV2 

1,2 


Ne IV] A2423 

2422.56 



0 

41278.89 

UV1 

5 


Ne IV] A2425 

2425.14 



0 

41234.43 

UV1 

5 


O II] A2471.0 

2470.97 



0 

40470.00 

1.01F 

3,6 


Oil A2471.1 

2471.09 



0 

40468.01 

1.01F 

3,6 

Mn IIA2577 

2576.88 

2.80E+08 

3.58E-01 

0 

38806.691 

UV1 

7,8,9 

Fe IIA2587 

2586.65 

8.94E+07 

7.17E-02 

0 

38660.054 

UV1 

1,2 

Mn IIA2594 

2594.50 

2.76E+08 

2.79E-01 

0 

38543.122 

UV1 

7,8,9 

Fe IIA2600 

2600.17 

2.35E+08 

2.39E-01 

0 

38458.993 

UV1 

1,2 

Mn IIA2606 

2606.46 

2.69E+08 

1.96E-01 

0 

38366.232 

UV1 

7,8,9 

Fe II 

* A2613 

2612.65 

1.20E+08 

1.22E-01 

384.7872 

38660.054 

UV1 

1,2 

Fe II 

* A2626 

2626.45 

3.52E+07 

4.55E-02 

384.7872 

38458.993 

UV1 

1,2 

Fe II 

* A2632 

2632.11 

6.29E+07 

8.70E-02 

667.6829 

38660.054 

UV1 

1,2 

Mg IIA2796 

2796.35 

2.60E+08 

6.08E-01 

0 

35760.88 

uvi 

10,11 

Mg IIA2804 

2803.53 

2.57E+08 

3.03E-01 

0 

35669.31 

uvi 

10,11 

Mg IA2853 

2852.96 

4.91E+08 

1.80E-00 

0 

35051.264 

UVI 

10,11 

Ti IIA3067 

3067.24 

3.47E+07 

4.89E-02 

0 

32602.626 

5 

10 

Ti IIA3074 

3073.86 

1.71E+08 

1.21E-01 

0 

32532.354 

5 

10 

He IA3189 

3188.67 

5.64E+06 

2.58E-02 

159855.97 

191217.05 

3 

12 

Ti IIA3230 

3230.12 

2.93E+07 

6.87E-02 

0 

30958.585 

2 

10 

Ti IIA3243 

3242.92 

1.47E+08 

2.32E-01 

0 

30836.425 

2 

10 

Ti IIA3385 

3384.73 

1.39E+08 

3.58E-01 

0 

29544.454 

1 

10 

[O II] A3727 

3727.10 



0 

26830.57 

IF 

3,6 

[O II] A3730 

3729.86 



0 

26810.55 

IF 

3,6 

[Ne III] A3870 

3869.86 



0 

25840.72 

IF 

13 

He I A3890 

3889.74 

9.47E+06 

6.45E-02 

159855.97 

185564.60 

2 

12 

[Ne III] A3969 

3968.59 



642.876 

25840.72 

IF 

13 

Ca IIA3935 

3934.77 

1.40E+08 

6.49E-01 

0 

25414.40 

1 

14 

Ca IIA3970 

3969.59 

1.36E+08 

3.21E-01 

0 

25191.51 

1 

14 

[S II] A4070 

4069.75 



0 

24571.54 

IF 

15 

[O III] A4364 

4364.44 



20273.27 

43185.74 

2F 

3,6 

He I A4473 

4472.76 

2.46E+07 

1.23E-01 

169086.91 

191444.48 

14 

12 


O III] A4960 

4960.30 



113.178 

20273.27 

IF 

3,6 


O III] A5008 

5008.24 



306.174 

20273.27 

IF 

3,6 


N I] A5199 

5199.35 



0 

19233.177 

IF 

3,4 


N I A5202 

5201.70 



0 

19224.464 

IF 

3,4 

Mg IA5169 

5168.76 

1.13E+07 

1.35E-01 

21850.405 

41197.403 

2 

10,11 

Mg IA5174 

5174.12 

3.37E+07 

1.35E-01 

21870.464 

41197.403 

2 

10,11 

Mg IA5185 

5185.05 

5.61E+07 

1.36E-01 

21911.178 

41197.403 

2 

10,11 

He IA5877 

5877.29 

7.07E+07 

6.10E-01 

169086.91 

186101.55 

11 

12 

Na IA5892 

5891.58 

6.16E+07 

6.41E-01 

0 

16973.366 

1 

11 

Na IA5898 

5897.56 

6.14E+07 

3.20E-01 

0 

16956.170 

1 

11 


O I] A6302 

6302.05 



0 

15867.862 

IF 

3,6 


O I] A6366 

6365.54 



158.265 

15867.862 

IF 

3,6 


O I] A6394 

6393.50 



226.977 

15867.862 

IF 

3,6 


N II] A6550 

6549.86 



48.7 

15316.2 

IF 

3,4 


N II] A6585 

6585.27 



130.8 

15316.2 

IF 

3,4 


He I A6680 

6680.00 

6.37E+07 

7.10E-01 

171134.90 

186104.96 

46 

12 


S II] A6718 

6718.29 



0 

14884.73 

2F 

15 


S II] A6733 

6732.67 



0 

14852.94 

2F 

15 


He IA7067 

7067.20 

2.78E+07 

6.95E-02 

169086.91 

183236.79 

10 

12 


Ar III] A7138 

7137.77 



0 

14010.00 

IF 

16 


O II 

A7322 

7322.01 



26810.55 

40468.01 

2F 

3,4 


O II 

A7332 

7331.69 



26830.57 

40470.00 

2F 

3,4 


O II 

A7333 

7332.75 



26830.57 

40468.01 

2F 

3,4 


a References: 1. IFuhr fc W iese (2006), 2. Nave & Johansson (2013), 3. I Galla gher & Moore d 1993D , 
4. iWiese Sz Fuh r (2007), 5. Kramida et al. (1999), 6. Wiese et ah| (I1996D , 7. IKr amida & Sansonetti 
(120131) , 8. [Aldcnius (2009), 9. [Den Hartog et al. ( 20 TU) , 10. iRuf foni & Pickering (2010), 
11. Kcllchcr & Podobcdova (2008), 12. Wiese & Fuhr (2009), 13. Kramida & Nave (2006), 14. 
I Safron ova & Safronova (2 01 ID , 15. IPodobe dova et al. ( 20091), 16. ISalomanl (120101) . 
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etc. are indicated by s, p, d, f, g, h, i and k, and the parentheses are often omitted. For example, the electron 
configuration of the ground state of Fe + is given by ls 2 2s 2 2p 6 3s 2 3p 6 3d 6 4s. In practice, we only show the last 
two or three subshells, e.g., 3d 6 4s for Fe + . 

o Spin-orbit ( L-S ) coupling: In L-S coupling, a given combination of (L, S) defines an atomic term. A given set 
of (L, S, J) defines an atomic level, and is indicated by 

2S+1 L P j, (A2) 

where p represents parity and is blank for even parity and “o” for odd. I Cohen et al.; ( 1987 ) defines an atomic 
state as a level with a given set of L, S, J and Mj. We find this definition is not used frequently in astronomy. In 
this paper, we use term and state interchangeably. We suggest to use microstate or quantum state for a given set 
of L, S, J and Mj (or L , S, Ms and Ml). With our choices, the ground term (state) of Fe + has L = 2 (D) and 
S = 5/2 with configuration 3d 6 4s and even parity, and is split into five levels with J = 1/2, 3/2, 5/2, 7/2, 9/2, 
among which the lowest (ground) level is 6 D g / 2 . The first excited term (state) of Fe + (that is relevant in the 
NUV spectroscopy) has L = 2 (D) and S = 5/2 with configuration 3d 6 4p and odd parity, and is split into five 
levels with J = 1/2, 3/2, 5/2, 7/2, 9/2, among which the lowest level is 6 Dg , 2 . Note each level has a degeneracy 
of g = 2J + 1 and the total multiplicity of a given term is (2 S + 1)(2J + 1). 

For a multiple-(sub)shell system, there are multiple terms with the same L and S for the outmost electron(s) 
due to different terms of electrons in the inner (sub)shells. iMoorel (1194511 labeled different terms with a front 
lowercase letter “a” - “z” to distinguish them and this labeling method has been used since. We also use this 
convention, with the labeling letters extracted from the NIST-ASD database. 

o Resonance: We refer to a transition between two different terms as a resonant transition if and only if the lower 
level of the transition is the lowest level of the ground term. For a permitted n on-resonant tra nsition, we label 
it with a right superscript asterisk This definition is different from that bv [Morton] (12003H . who considered 
transitions whose lower term is the ground term as resonant transitions, regardless of which level. Note we only 
apply this to permitted transitions, so forbidden (semi-forbidden) lines like [O III] A5008 (C II1A2326), although 
the lower energy level is not the lowest level, are labeled with the brackets “[]” (“]”) onl^Q. We refer to a 
transition between two different levels, i.e., with different J-values, within the same term as a fine-structure 
transition. Fine-structure transitions are forbidden since the configuration does not change and thus A l = 0. 

o Multiplet groups: IMoorel (1945) labeled the multiplets first in the order of increasing excitation energy of the 
lower term then the higher term. At wavelength A > 3000 A, the multiplets were inclu ded in the Revise d 
Multiplet Table ( Moore!fl972 ). and at shorter wavelength in the Ultraviolet Multiplet Table (Mooreil 1950 1 11952 1). 
Though the data have improved over the years, this convention has been used by many author s. Newly-identified 
multiplet groups were usually added between the original groups by a decimal system (e.g., iGallagher fe Moorel 
119931) and forbidden lines are denoted with “F” following the Multiplet Number. We notice that, in practice, the 
Multiplet Numbers are rarely used for transitions in the optical. Nonetheless, for consistency with th e lit erature, 
we follow the same convention, and present the Multiplet Numbers taken from IMoorel ( 1950 . 19521 fl972t ) or new 
references we found in the literature. 


• Ionization potentials 

Table ED presents the ionization potential for elements with atomic number up to 30, in unit of eV. We have 
included the most abundant isotopes ionized up to seven times. The columns N p and N n are the numbers of protons 
and neutrons. When available, we show seven significant figures. 

• Line list 

In Table (A2l we list the lines associated with star-forming galaxies between 2200 A and 7500 A. We detect most of 
them in the composite spectra of the emission-line galaxies from the eBOSS pilot observations. We do not list the 
categories of origin and refer the reader to Section 13.31 for a summary of typical origins for the lines in the NUV. 
We present vacuum wavelength and some of the lines are therefore labeled as 1A longer than if labeled with the air 
wavelength, such as [O III] A5008. For useful formulae for the conversion between vacuum and air wavelengths, we 
re fer the reader tolCiddorl (jl996l) . We have omitted the hydrogen lines, a recent compilation of which can be found 
in IWiese fc Fuhr] (2009f), and the G4300-band, which is a blend of many lines predominantly by CH and Fe I (e.g., 
IWorthev et al.lll994H. 

For permitted transitions, we list the Einstein A coefficient (A u i) and oscillator strength (fi u ), which are essential 
quantities for the radiative transfer analysis presented in this paper. For all the transitions, we list the energies of 
the lower (E{) and higher levels (E u ), relative to the ground level (of the ground term). The energies are in unit of 


26 According to the selection rules, the transitions in Fe II UV4 
and UV5 shown here, A2250, A2261, (*)A2270 and (*)A2281 are 
semi-forbidden. Because the selection rules apply to emission lines 
and historically Fe II lines were more often observed in absorption, 


these lines are usually not marked with “]”. Nonetheless, their 
transition probabilities are several orders-of-magnitude higher than 
other semi-forbidden lines. We treat them as permitted transitions 
here. 
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Fig. A1.— The energy-level diagrams for the Fe II UV2 and UV3 groups. The numbers next to the arrows represent the vacuum 
wavelength, the Einstein A coefficient, and the oscillator strength. 

cm -1 , and can be converted to eV when multiplied by 1.2398E-4. The transition wavelength (in A) can be derived 
with 1E8 /{E u — Ei). 

• Energy-level diagrams 

The relationships between the different resonant and non-resonant channels are essential to our analysis of radiative 
transfer processes. We illustrate these relationships with the energy-level diagrams of the corresponding multiplet 
groups. The diagra m fo r the Fe II UV1 group is presented in Sec tion HOI (Figure 0. Here we present those for Fe II 
UV2-UV5 fFigures lAll fc lA2l) . Mg II UV1 and Mg I UV1 (Figure lA3]h We also include C II UV0.01 and Mn II UV1 
(Figure |A4|, and Ti II No. 1, 2 & 5 f Figures I A5ll as they merit further attention once we obtain more data in the near 
future. For other lines identified in the composite spectra of ELGs, such as [O II] AA2471.0,2471.1, their energy-level 
diagrams can be found in standard textbooks (e.g., Appendix E in [Drain ell201 111 . 

The terms and symbols were introduced in the nomenclature above. For each transition, we show the vacuum 
wavelength, the Einstein A coefficient, and the oscillator strength. The black lines with upward arrows represent 
the observed absorption lines. The green lines with downward arrows show the observed or expected emission lines, 
with solid indicating permitted lines, dashed semi-forbidden (C II]). The brown lines with downward arrows represent 
fine-structure transitions. 

B. COMPARISON OF ELGS WITH LOCAL STAR-FORMING REGIONS 

In Section 13751 we compare the ELG composite spectrum (in blue) with the composite of the nine local SF regions 
(in red) at 2200 A < A < 2900 A in Figure 0 In Figure IB11 we zoom in on the most prominent features to show 
the details, as we did in Figure 0 for the comparison with quasar absorption-line systems. Because the UV atlas 
by iLeitherer et al.l (|2011f ) does not have wavelength coverage longer than 3200 A, we do not compare the spectra at 
3000 A < A < 4000 A as i n the botto m panel of Figure 0 As mentioned in Section 1331 we have corrected the original 
wavelength in ILeitherer et all (1 20 1 ill by —0.7 A based on the position of C III A2298. The main differences are in the 
non-resonant emission lines and the absorption line ratios. See the main text for details. 

C. THE DEPENDENCES OF THE NUV FEATURES ON [O II] EQUIVALENT WIDTH AND LUMINOSITY 

The wavelength coverage of eBOSS allows us to study the dependences of the NUV features in the ELG spectra on 
the [O II] AA3727,3730 properties. We divide the NUV sample of 8620 ELGs into half subsamples based on the total 
rest equivalent width and luminosity of the [O II] AA3727, 3730 doublet. In Section 0 we discuss in depth the emission 
strength, the observed absorption profiles of Mg II, the emission-corrected absorption strength, and the unified velocity 
profiles. He re w e present some additional details of the analysis. 

In Figure ED we present the distributions of total rest equivalent width (IUjq 3 /^ 8 ) and luminosity (in logarithmic 

scale, log 10 Ljqjj 8 ). The median value of is about 51.4 A and the median of log 10 Lpij 8 is about 41.6dex. 
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Fig. A2.— The energy-level diagrams for the Fe II UV4 and UV5 groups. 
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Fig. A3.— The energy-level diagrams for the Mg II UV1 and Mg I UV1 groups. 

Figure IC2l presents the observed velocity profiles of the non-resonant Fe II* emission. The dependence on the [O II] 
rest equivalent width is larger than on the luminosity. 

We show examples of the observed absorption profiles for the different subsamples in Figure lC3l The observed profiles 
are affected by the scattered emission filling in on top of the true absorption. The largest correlation is between the 
Mg II pr ofile s and the [O II] rest equivalent width, which we discuss in Section [5] 

Figure IC4l pres ents the emission-corrected absorption profiles, based on the observation-driven iterative method de¬ 
scribed in Section l4~2l We show the correlations between the rest equivalent width of the emission-corrected absorption 
with the [O II] AA3727, 3730 properties in Figure [HI 
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Fig. A4.— The energy-level diagrams for the Mn II UV1 and C II UV0.01 groups. Note although we show the Einstein A coefficient and 
the Oscillator strength of the C II transitions, they are very low for these semi-forbidden lines. 
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Fig. A5.— The energy-level diagrams for the Ti II No. 1, 2 & 5 groups. 
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Fig. Bl.- The composite continuum-normalized spectrum of 8620 ELGs at 0.6 < z < 1.2 (blue) in eBOSS, compared with the composite 
spectrum of local star-forming regions (red) taken with the FOS/GHRS on HST. 
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Fig. Cl.— The distributions of the total [O II] AA3727, 3730 rest equivalent width (left) and luminosity (right, in logarithmic scale). To 
calculate the luminosity, we assume the ACDM cosmological model. 




Fig. C2.- The dependences of the emission velocity profiles on the [O II] AA3727, 3730 rest equivalent width (left) and luminosity (right). 
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Fig. C3.— The dependences of the observed absorption velocity profiles on the [O II] AA3727, 3730 rest equivalent width (left) and 
luminosity (right). 
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Fig. C4.— The dependences of the emission-corrected absorptio n ve locity profiles on the [O II] AA3727 , 37 30 r est equivalent width (left) 
and luminosity (right). The color scales are the same as in Figure [To] based on the orders given by Eas. 113l and Il4l 




























